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Some of the outstanding questions of particle physics today concern the neutrino sector, in par¬ 
ticular whether there are more neutrinos than those already known and whether they are Dirac or 
Majorana particles. There are different ways to explore these issues. In this article we describe 
neutrino-mediated decays of charged pseudoscalar mesons such as tC, K ± and B ± , in scenarios 
where extra neutrinos are heavy and can be on their mass shell. We discuss semileptonic and lep- 
tonic decays of such kinds. We investigate possible ways of using these decays in order to distinguish 
between the Dirac and Majorana character of neutrinos. Further, we argue that there are significant 
possibilities of detecting CP violation in such decays when there are at least two almost degenerate 
Majorana neutrinos involved. This latter type of scenario fits well into the known neutrino min¬ 
imal standard model (j/MSM) which could simultaneously explain the Dark Matter and Baryon 
Asymmetry of the Universe. 

PACS numbers: 14.60St, 11.30Er, 13.20Cz 


I. INTRODUCTION 

To date it is unclear whether the neutrinos we know are Dirac of Majorana fermions. Unlike Dirac fermions, 
Majorana fermions cannot be distinguished from their own antiparticles. As a consequence, processes involving Dirac 
neutrinos conserve charges such as Lepton Number, while processes involving Majorana neutrinos will not conserve 
them. There exist several processes which may clarify the Majorana or Dirac nature of neutrinos. Among such 
processes the most prominent are neutrinoless double beta decays (0isf3/3) in nuclei [TUTU]. Other such processes are 
specific scattering processes [TTMT] . and rare meson decays [22U40] . 

A related issue in neutrinos physics is the absolute mass values of the known neutrinos. While the experimental 
evidence of neutrino oscillations within the known three flavor states [HfHSl clearly shows that these particles cannot 
all be massless, the oscillations are only sensitive to mass differences, not to their absolute values. In contrast, Ou/3/3 
decays are sensitive to the absolute mass and may help in their determination, if neutrinos turn out to be Majorana 
particles. So far the best bounds on the absolute masses of the light neutrinos come from Cosmology m„ > 0.23 eV 

S3- 

A pending question is then why light neutrinos are so light, specifically so much lighter than all other Standard 
Model (SM) fermions. Interesting enough, the existence of such very light neutrinos can be explained via the seesaw 
mechanism |48H52| where more neutrinos are required and where all of them are, in general, Majorana particles. In 
the simplest form of this mechanism, the masses of the light neutrinos are ~ M 2 d /Mr (<1 eV), where M D is an 
electroweak scale or lower. At the same time, additional neutrinos, usually much heavier (with masses Mr 3> 1 
TeV) and sterile under electroweak interactions except through small mixing with the SM flavors, are required. This 
mixing is suppressed as ~ Mr/Mr (<Cl). Besides the simplest scenario, there are other seesaw scenarios in which 
the heavy neutrinos may have lower masses, namely near or below 1 TeV |53h59j and even near the 1 GeV scale or 
below [T 51 filO'Wihj , and at the same time their mixing with the SM flavors may not be so extremely suppressed as in 
the original scenarios. 

In the first part of this work we discuss lepton number violating (LNV) semileptonic decays of charged pseudoscalar 
mesons such as K ± and B ± , mediated by a heavy Majorana neutrino on its mass shell, cf. Ref. [261. The pions, which 
are the lightest mesons, can have only leptonic decays; we discuss hypothetical leptonic decays that could be mediated 
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by on-shcll heavy neutrinos. Such decays could be either lepton number conserving (LNC) or lepton number violating 
(LNV) when the mediating neutrinos are Majorana particles, cf. Ref. [37], while only LNC decays occur when the 
neutrino is of Dirac type. We present ways of determining the nature of neutrinos using the differential decay rates 
of these processes. 

Yet another interesting issue in neutrino physics is the possible existence of CP violation in the lepton sector, 
a phenomenon that could be measured, for example, in neutrino oscillations [66]. Alternatively, as we present in 
the second part of this work, leptonic CP violation may show in leptonic LNC and LNV decays of charged pions, 
cf. Ref. [32], or in semileptonic LNV decays of and R , cf. Refs. [3HSS]- It turns out that such CP violation 
becomes appreciable and possibly detectable in these decays if the scenario contains at least two on-shcll heavy 
neutrinos that are almost degenerate. Interestingly, this scenario fits well into the neutrino minimal standard model 
(z/MSM) L60, SS[S7H72I which contains two almost degenerate Majorana neutrinos of mass ~ 1 GeV and another lighter 
neutrino of mass ~ 10 1 keV, besides the three light neutrinos of mass <1 eV. This model can explain simultaneously 
the existence of neutrino oscillations, dark matter and baryon asymmetry of the Universe. Furthermore, in more 
general frameworks of low-scale seesaw, baryon asymmetry (but not dark matter) is explained while keeping even 
larger values of the heavy-light mixing m than in the i/MSM, and in such frameworks the case of almost degenerate 
Majorana neutrinos is preferred m since it allows larger mixings. CP violation effects in the neutrino sector in 
scenarios with nearly degenerate heavy neutrinos have also been investigated earlier [75[ 76] using a more detailed 
formalism, although it amounts to the same effect described here which is the interference of amplitudes with two 
slightly different dispersive and absorptive parts in the neutrino self energy. 

In Section [n] we discuss the LNV semileptonic decays of mesons M ± —> £f£fM lzf , mediated 
by an on-shell Majorana neutrino (which henceforth we call N) where M is a heavy pseudoscalar 
(M = K, D, D s , B, B c ), M is a lighter pseudoscalar, and £.j (j = 1,2) are charged leptons, and we present there the 
corresponding branching ratios. In Section [TTI| we present the expressions and values of the branching ratios for the LNC 
and LNV leptonic decays of charged pions mediated by on-shell N sterile neutrino, 7r =l= —► e ± N —> e^e^fjFv, as well as 
the differential branching ratio dQr/dE^ for these decays. We discuss the possibilities of detecting such branching ratios 
and to discern from them the Majorana or Dirac nature of neutrinos. In Section [TV] we then extend the analysis of the 
mentioned leptonic and semileptonic decays to a scenario where we have at least two heavy on-sliell sterile neutrinos in¬ 
volved (N\, N 2 ), and we present an analysis of CP-violating asymmetries Acp = [T{M~) — r(M + )]/[r(M“) +r(M+)] 
for such processes. In Appendices 1 and 2 we present explicit formulas for our LNV semileptonic decays, and in Ap¬ 
pendices 4 and 5 explicit formulas needed for the analysis of our LNC and LNV leptonic decays of the charged pion. 
Appendix 3 contains formulas needed for evaluation of the decay width of the heavy neutrino N, and in Appendix 6 
we derive an identity relevant for CP violation asymmetry. In Section [V] we discuss and summarize our results. 


II. LEPTON NUMBER VIOLATING SEMILEPTONIC DECAYS OF SCALAR MESONS 

If there is a Majorana sterile neutrino N with mass M jv ~ 1 GeV, its existence could be discerned by detecting 
semileptonic LNV decays of heavy mesons mediated by on-shell N. Here we will consider such LNV decays M± —>■ 
£f N —> where M and M are pseudoscalar mesons (M = K, D, D Sl B, B c ; M = n, K, D, D s ) while £\ and 

i 2 are charged leptons (e, // or r), cf. Figure [IJ 

Large part of this Section uses the results of Refs. [21133, and for certain general formulas, those of Refs. [31 EDI- 

We consider a scenario where we have at least one heavy sterile neutrino N , which has (small) mixing B^ with 
the three known neutrino flavors ve {£ = e, fx, t) 
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Vi = ^2 Bivk Vk + Bfi N N + ... ( 1 ) 

k =1 

Here, (k = 1,2,3) are the three light mass neutrino eingenstates. The considered decays may be appreciable only 
if N can go on its mass shell (in the s-type channel, Figure [l]), creating a very large resonant enhancement of order 
totv/Pjv, condition which is fulfilled if 

(M m > + M t2 ) < M n < {M m - M tl ) , or/and 

{Mm 1 + M ( t ) < M n < {M m ~ M ^ 2 ) (2) 

Consequently, only tree level resonant amplitudes need to be considered. The mixing matrix B ap¬ 
pearing in Equation |l]) should be unitary, implying that the PMNS 3x3 block Bt Vk {£ = e,n,T 

and k = 1,2,3) is in general not unitary. If one adds extra and heavy neutrinos—as in most seesaw 
models—the unitarity of B thus provides upper limit constraints on the heavy-to-light mixing elements EH IZ23-ES3- 
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This is in part one of the reasons for the high suppression suffered by all lepton flavour violating processes involving 
heavy neutrinos. 




FIG. 1: The lepton number violating (LNV) semileptonic decay M + (pm) —> M ( pm ') mediated by a Majorana 

neutrino N: (a) the direct (D) channel; (b) the crossed ( C ) channel. 


A. Branching Ratio for M 3 


IflfM T 


The decay width of the considered decays can be written as final particles’ phase space integral of the square of the 
reduced decay amplitude T{M ± ) (summed over helicities of charged leptons) 


T{M ± = {2-5 tlh ) 


1 1 


1 


±\|2 


2! 2 M m (2?r) 5 


d 3 |T(M ± )| 


(3) 


Factor 1/2! above is the symmetry factor when the two produced leptons are equal; d 3 is the integration differential 
of the final three-particle phase space 


d 3 = 


d 3 pi d 3 p 2 


d 3 PM 


; —c£ (4) (Pm ~ Pi ~ P2 ~ Pm ') 


2-E^i (pi) 2 Ee 2 (p 2 ) 2Sm' (pm> ) 

The resulting decay width can be written as 

T(M ± —>• ef = (2-Se ie2 )\k\ 2 \f{DD*) + f{CC*)+f ± (DC*)+f ± {CD*) 

where k is the mixing factor 


(4) 


(5) 


k = 


( 6 ) 


and r±(Xy*) are the normalized (be., without the explicit mixing) contributions from the A' channel and the 
complex-conjugate of the Y channel (. X,Y = D,C, where D and C stand for the direct and crossed channels) 


T±(Ay*) = K 2 I J d 3 P(X)P(Y)*M 2 n T ± (X)T ± (Y)* 


(7) 


The expressions for T^(A)T^(y)* ( X,Y = D,C ) are given in Appendix 1. T±(X) is the relevant part of the 
amplitude in the X channel and forms part of the total decay amplitude T(M ± ), cf. Appendix 1. In Equation (J5|, 
notice that subscripts ± for the contributions T(DD*) and T(CC*) are unnecessary because \T + (D)\ 2 = |T_(D)| 2 
and |T + (C)| 2 = |T_(C)| 2 . P(X) (X = D, C) are the propagator functions of the intermediate neutrino N in the two 
channels 


[(pm -pi) 2 - Mh +iT N M N ) 
[(Pm -P 2) 2 ~ Mh +iT N M N ) 

The overall constant K 2 in Equation Q is 


(8a) 

(8b) 


K 2 = G 4 F f 2 M f 2 M'\V QuQd V quqd \ 2 


( 9 ) 
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Here, /m and /m' are the decay constants of M ± and M U and VQ u Q d and U 9tigd are the corresponding CKM matrix 
elements. We denote the valence quark content of M + as Q u Qd\ of M + as q u qd- 

When the intermediate neutrino N has such a mass that it is on mass shell. Equation (J 2 J , the squares of the 
propagators ((81) are reduced to Dirac delta functions because < Mjv 


\P(X)\ 2 = 


( Pm ~Pk ) 2 - Mjv + «r at Mat 

5{{Pm ~ Pk) 2 — Mh) (T n <^M n ) 


m n t n 


( 10 ) 


where pk = Pi,P2 for X = D,C. In this on-shell case, the DD* and CC* contributions in Equation ([5]) are large, 
and the interference contributions DC* and CD* are negligible in comparison (c/. Ref. [39J for details on this point), 
leading to 


r(M ± -»■ = (2 - 6e li2 )\k\ 2 T{DD*)+T(CC*) 

= r(M ± -> ef£%M'*-,DD*) +T(M ± 


£f£fM T -,CC*) 


(11a) 

(lib) 


when £ x = £ 2 , we even have E (DD*) = T(CC*). The normalized decay width T(DD*) can be calculated explicitly, 
and it turns out to be 


f (DD-) = if 

1287T 2 r N \ VN VN J 

and T(CC*) is obtained by the simple exchange yg 1 O yz 2 

r(CC*) = f (DD*)(y £l ^ y, 2 ) 


QiyNWt^yt^y') 


The notations used in Equations ([12]) and (131 are 

2/ 2j 2 / 3 ) = y\+y\ + y\-2y x y 2 -2y 2 y- i ~2y z yi 

VN = > Vis = TFT > V' = T75T . (4» = ti,h) 


M 


M 


Mm 


M 2 m 


( 12 ) 

(13) 

(14a) 

(14b) 


and the function Q(yN',ye 1 ,ye 2 ,y r ) is given in Appendix 2. In the limit of massless charged leptons (y^ = y ^ 2 = 0), 
the expression |l2|) reduces to 


T(DD*)\ 


K 2 M 5 m M n 


M tl —M t2 = 0 


256tt 2 T n 


Vn 


(15) 


We note that the expression ( |l2| , although having the explicit mixing dependence factored out [cf. Equation ©]. 
contains the dependence on the mixing coefficients Bpjy in the denominator d ue t o the N-decay width there T^v oc 
\B(n\ 2 {(■ = e,/z,r, see below). This factor 1 /Tjv in T(DD*) of Equation (12) represents the 7V-on-shell effect 
Equation (10). As a result, the considered width —> £ x £ 2 M T ) is by many orders of magnitude larger when 

N is on-shell than it would be if N were off-shell. For more quantitative analyses, it is thus important to have an 
expression for Tjv as a function of mass Mn- Using the results of Ref. [39], we can write this decay width as 


T/v = KY n (M n ) 

where the corresponding canonical (i.e., without any mixing dependence) decay width is 

GlM% 


Tjv(Mjv) = 


967T 3 


(16) 


(17) 


and the factor K, contains the dependence on the heavy-light mixing factors 

/C(Mjv) = K = Af e N |-B e jv| 2 +Af(iN \B^n\ 2 +A f T N | B t n 1 2 (18) 

In this expression, J\f^ N {M N ) = Afg N (£ = e,/q r) are the effective mixing coefficients; these are numbers ~ 10°-10 1 
which depend on the mass Mn- In Appendix 3 we write down the relevant formulas for the evaluation of these 
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coefficients. The results of these evaluations are presented in Figure [2j for the case of Majorana and Dirac neutrino 
N, in the entire neutrino mass interval 0.1 GeV < Mm < 6.3 GeV which will be of interest in this work. For further 
clarifying remarks we refer to Appendix 3. Equations (11) and ( Jl2| imply that T(M ± —> CfifM ff) is proportional 
to 1//C (oc l/\B m \ 2 ). Hence, we can define a canonical branching ratio Br, being the part of the branching ratio 
Br(DD*) = T(M ± —> =F ; DD*)/T(M ± —► all) with no explicit or implicit heavy-light mixing factors 


Br(M ± ->• efefM T -DD*) = 


T(M ± £?££M T ;DD*) _ (2 - S tl 


= ( 2 -, 


r(M± all) 

-2Br(2/iv; y tl ,yt 2 \ y') 


T(DD*) 


F(M± ->■ all) 

W,) T 

Use of the expressions (12) and Equations (16) and then gives for the canonical branching ratio the following 


(19a) 

(19b) 


expression: 


Br (DD*) = Bv(y N -y ei ,ye 2 ,y') 
3tt I< 2 M m 


8 G|r(M± -> all) y% 


\-X 1/2 (l,y N ,y ei )X 1/2 ( 


-i V V(-2 \ f\ 

, i, —. — Q{yN m ,yii,ytz,y ) 
V VN VN 


( 20 ) 


where the notations (141 and © are used. In the limit of massless charged leptons (M^ = M( 2 = 0) this expression 
becomes simpler 


Br(yAr; 0, 0, y') = 


37r 


K 2 M m 


16 G 2 f T(M± 


all) 


(1-Wv) 2 1- — 
VN 


( 21 ) 




FIG. 2: The effective mixing coefficients Min {£ = e,y,,r) appearing in Equations (161 (181, as a function of the mass Mn of 
the neutrino N. The left-hand figure is for the case of Majorana neutrino; the right-hand figure for the case of Dirac neutrino. 
For details we refer to Appendix 3. 


B. Effect of the Long Neutrino Lifetime on the Observability of —> £^£^M T 


In the mentioned branching ratios, an often important effect of suppression due to the decay (*.e., nonsurvival) 
probability was not included. Namely, if the detector for the considered decays has a certain length L, the produced 
(on-shell) massive neutrino N could survive during its flight through the detector, and would decay later outside it. 
Such decays are thus not detected and should be eliminated from the width and the branching ratio of the considered 
process M ± —»• ifd^M ff, by introducing a suppression factor (nonsurvival probability) Py = 1 — exp[— £/(tjv7jv)]> 
where t ~ L//3m is the time of flight of N through the detector (/3jy is the velocity of N in the lab frame), and 
Tat = (1 — is the Lorentz time dilation factor. Hence, the suppression factor, which should multiply the 

branching ratio, is 


P N = 1 — exp 


L 

TnInPn 


1 — exp 


LT n 

7iv 


( 22 ) 
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In the last relation, we used /3n ~ 1 and tjv = 1/Tjv [= l/r(IV —► all)], in the units used here (c = 1 = h). This 
decay-witliin-the-detector probability Pn has been discussed and presented for the processes with intermediate on- 
shell particle (such as N) in Refs. 15 } [57H551 I80H812] . In this respect, here we follow mostly the notations of Ref. [391 . 
Usually, the quantity P/v is small and is then written as 


Pn ~ L/(tn^nPn) {~L/{t n 7 jv)) if P/v <C 1 

which agrees with Equation (22) in the limit of small Pn- The suppression factor (|22| can be rewritten as 


■ N 


= 1 — exp ( — 


L 

Ln 


= 1 — exp ( — /C 
Ln 



Ln 


if P N < 1 


(23) 


(24a) 

(24b) 


Here, is the decay length, and is the canonical decay length (independent of the mixing parameters N ) 


J N 

--1 


= LnK 


l n — 


g 2 f m 5 n 

967t 3 


(25a) 

(25b) 


T n {M n ) = 

7 n 7 n 

where K, and Tn(Mn) are from Equations (16)—(18), cf. also Figure^ Equation (24b) suggests that it is convenient 
to define a canonical (i.e., independent of mixing) probability Pn for the decay of N within the detector as 


— lm — / L \ ~ 

N = In ^ P#KP Hlm ) (26) 

We present the inverse canonical decay length, L ^, for 7 n = 2, in Figure[i]as a function of Mn- We note that Ln 
increases very fast (as M FT ) when Mn increases. Therefore, the supression due to the factor Pn may not necessarily 
be strong (i.e., Pn it. 1) for semileptonic LNV decays of heavier mesons M , such as B ± . If we use eyeball estimates 



Mjv(GeV) 




FIG. 3: The inverse canonical decay length L N * = rjv(Apv)/7.iv, Equation (25b I, in units of inverse meters (m x ), as a function 
of the neutrino mass Mn, with the Lorentz lab time dilation factor chosen to be yjv 1= (1 — = 2. The y axis can also 

be interpreted as the canonical probability Pn = (lm)/Ljv (dimensionless), Equation (26). 


for the coefficients A r m of the left-hand Figure [2J approximate expressions for the factor K, of Equation ( 18 ) for 
Majorana neutrinos can be written 


K s 

s 15 \B eNj 

\ 2 + 8\B^ Nj \ 2 + 2\B tNj \ 2 

(K decays) 

(27a) 

K. s 

« 7(1 B eNj \ 

2 + \BnNj 2 ) + 2 \B TNj 2 

(. D , D s decays) 

(27b) 

K. p 

« 8(1 B eNj \ 

\ 2 + \B, N] \ 2 ) + 3\B tNj \ 2 

(. B , B c decays) 

(27c) 
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TABLE I: Presently known upper bounds for the squares |-£W | 2 of the heavy-light mixing matrix elements, for various specific 
values of Mjv. 


Mjv [GeV] 

IBeJVr 

\B,n\ Z 

| BrirT 

0.1 

(1.5 ±0.5) x 10~ s 

(6.0 ±0.5) x 10~ b 

(8.0 ±0.5) x 10 ' 4 

0.3 

(2.5 ±0.5) x 10~ y 

(3.0 ±0.5) x 10~ M 

(1.5 ±0.5) x 10 i 

0.5 

(2.0 ±0.5) x 10~ M 

(6.5 ±0.5) x 10~ y 

(2.5 ±0.5) x 10 ~ 2 

0.7 

(3.5 ±0.5) x 10~ s 

(2.5 ±0.5) x 10~ Y 

(9.0 ±0.5) x 10~ a 

1.0 

(4.5 ±0.5) x 10~ s 

(1.5 ±0.5) x 10~ Y 

(3.0 ±0.5) x 10 ;i 

2.0 

(1.0 ±0.5) x 10 ~ 7 

(2.5 ±0.5) x 10 ~ 5 

(3.0 ±0.5) x 10 ~ 4 

3.0 

(1.5 ±0.5) x 10~ Y 

(2.5 ±0.5) x 10~ b 

(4.5 ±0.5) x 10 ” 5 

4.0 

(2.5 ±0.5) x 10~ Y 

(1.5 ±0.5) x lO ' 13 

(1.5 ±0.5) x 10 “ 5 

5.0 

(3.0 ±0.5) x 10 ~ 7 

(1.5 ±0.5) x 10 ~ 5 

(1.5 ±0.5) x 10 -5 

6.0 

(3.5 ±0.5) x 10~ Y 

(1.5 ±0.5) x 10~ b 

(1.5 ±0.5) x lO ’ 5 


TABLE II: Rough estimates of upper bounds for |-B«jv | 2 {£ = e, /r, t), for Mjv in three different ranges around the values 0.25, 
1, 3 GeV; and the inverse of the canonical decay length, L N (in units of m~ 1 and for 7 jv = 2). 


Mjv [GeV] 

|ReJv | 2 

\B,n\ 2 

| B t n 1 2 

L n [m 1 ] 

« 0.25 

10 ~ s 

10 "' 

10" 4 

0.11 

« 1.0 

10~ 7 

10~ 7 

10 -2 

1.1 x 10 2 

« 3.0 

10~ 6 

10~ 4 

10" 4 

3 x 10 4 


In order to estimate better the values of K, Equations ([27]) and thus the suppression factor Pm Equation (|24j), we 
need to know the present upper bounds for the squares \B^m\ 2 as a function of M These upper bounds we take 
from compilation of values of Ref. [25] . based in turn on upper bound values obtained in Refs. [55H5S] . We present 
them in Table [TJ for specific chosen values of Mn in the mass range of interest. We remark that the upper bounds 
have in some cases strong dependence on the precise values of M n , see Ref. [29j for further details. In order to use 
only rough estimates for the values of V, we present in Table [il] order of magnitude values for upper bounds of \B^ N \ 2 . 
These rough upper bounds are given for three typical ranges of our interest: Mn around 0.25; 1; 3 GeV. They are 
relevant for the decays of I \; (D, D s ); (B, B c ), respectively. The corresponding values of the inverse of the canonical 

decay length, L^, are included. As seen in Tables |l| and |h| the upper bounds for |.B r jv | 2 are at present si gnific antly 


less stringent and are expected to become more stringent in the future. When we combine Equations (24b I with 
(27) and Table [nj we obtain for the decay-within-the-detector probability Pjv = P at/C the following estimates and 
upper bounds, relevant for the K decays {Mn « 0.25 GeV), D and D s decays (My « 1 GeV), and B and B c decays 
{Mn ~ 3 GeV), all when L = 1 nr and 7 n = 2: 


r(M jv « 0.25GeV) 

PS 

1.7] B eNj | 2 ± 0.9] B^ | 2 {+0.2\B rNj | 2 ) 



< 

10” 8 ±10" 7 (±10 -5 ) 

(28a) 

P n {M n « lGeV) 

PS 

0.8 • 10 3 | B eNj | 2 ± 0.8 • 10 3 | 2 (±2 • 10 2 | B tN . | 2 ) 



< 

10" 4 ± 10" 4 (±10°) 

(28b) 

P n {M n w 3GeV) 

RS 

3 • 10 5 |ReAr.] 2 ± 3 • lO 5 ]^] 2 (±1 • 10 5 |R t jv.| 2 ) 



< 

10° ± 10° (±10°) 

(28c) 


In order to have the analysis and the formulas simpler, in the rest of this Section we will assume that one mixing 
parameter, \Bin\, dominates over the other two mixing parameters: 


Bin I ^ \B, 


•e.'N\ 


{£ ± 1 ) 


(29) 


For example, it may well be that £ = /i, i.e., that \B^n\ \B e N\i \B t n\- Then, of the branching ratios Br(M ± -± 
the largest will be M ± —> £±£±M^ which, according to Equations (29) and (19) (note that DD* and 
CC* give the same contribution since i\ = £2 now), is: 


Br(M ± —»• UrUrM T ) =4 




(30) 


























TABLE III: Values of the factor Br e ff, Equation (321, with L = 1 m and 7 jv = 2, for some of the LNV decays M ± —> t' ± t' ± 7r :F . 
The value of Mn is chosen such that the maximal value of Br e ff is obtained (the value of Mn is given in parentheses, in GeV). 
For M ± = K ± , two different values are given, for £ = e and £ = fi. For all other cases, £ = fj, is taken (when £ = e the values 
are similar). 


M ± : 

'oT 

II 

u 

3. 

II 

■H 

D= 

Df 

B ± 

Bf 

Br e ff 

6.8 (0.38) 

3.8 (0.35) 

3.9 (1.39) 

70. (1.47) 

0.96 (3.9) 

199. (4.7) 


Multiplying this expression by the probability P/v of the decay in the detector, Equation (26), we obtain the effective 
branching ratio Br e ff 


Br eff (M ± -)• e^M'*) = P N Br(M ± = 

L 


Pn[ ^) IC 


4^^Br' 

K 


4|P^jv| 4 PjvBr = |P^jv| 4 ^ ^ Br ef f 


(31a) 

(31b) 


We see that in the effective branching ratio, Br e g, the complicated dependence on mixing parameters encoded in K. 
[cf. Equation (18)] disappeared because factors 1C cancel here. All the mixing effects in Br e ff are in the simple factor 
IBfArl 4 . Unfortunately, this factor represents a strong suppression, in comparison with Br of Equations (191 where 
Br oc \k\ 2 /K = \Bt^B l N \ 2 JK ~ \B m \*. 

In the last identity (31b) we introduced the canonical (i.e., without any mixing dependence) effective branching 
ratio Bi'eff 


Br ef r = 4PArBr = 4 ( Br 


(32) 


where we recall that Br was defined in Equations (19) and (|20|). Here Br e g is half the value of Br e g- in Ref. 


because the latter expression referred to the exchange of two Majorana neutrinos instead of one. 

Only when M ± = B ± or P;f, i.e., when the mass of the on-shell N can be high (Mn > 1 GeV), would it be possible 
to have P/v ~ 1 [Equation (28cl|; and in such a case Equations (31) do not apply, but rather Equations ( fTT(T| ) . i.e.. 
Br e ff = Br in this case. Figures [4j|7] show the effective canonical branching ratios ( |32] ) as a function of the neutrino 
mass Mn, for various considered LNV decays of the type M ± —» Figure [ 4 ] for M = AT; Figure [H^,, b 

for M = D,D S , respectively; Figures [ 6 ji and [7}i for M = B,B C , respectively. We took £ = e, //, and L = 1 m and 
7 n = 2. For the case when Pn ~ 1 and hence the estimates Equations (30) apply, Figures [ 6 ]d and[?]) present the 
branching ratios Br(:c) as a function of Mn, for M^ = B ± and B±, respectively. For the meson decay constants and 
CKM matrix elements, needed for the evaluation of I \ 2 factor of Equation (|9j), and for the masses and lifetimes of 
the mesons, we used the values of Ref. m- The values of the decay constants f b and fs c were taken from Ref. [98] : 
f b = 0.196 GeV, fs c = 0.322 GeV. We note that the presented formulas for Br eff and Br can be evaluated also for 
the decays M ± —> ifl^M t when £\ ^ £ 2 - Furthermore, when the final leptons are r leptons (and M ± = B ± or 
B±), the values of branching ratios turn out to be sim ilar to those in Figur es [ 6 ] and [?| but the range of Mn in this 

~~ displays values of Br e fj for representative values of Mn in 


Table 


III 


case is shorter: M m > + M T < Mn < Mm — M. 

t ^^.s < an^!liMratwc exan/ple', let us consider the decays Df 
proposed at CERN-SPS jSD., 811. and, in addition, let us assume that 




which is one of the preferred decay modes 
|P M zv| 2 is the dominant mixing. In such a case, 


Equations (31) and Table III imply for the experimentally measurable branching fraction Br 0 g 


Br e ff (Dj 


= P/vBr(P 4 


r 1 


) ~ 10 2 \B^n\ 


(33) 


The present rough upper bound on the mixing for Mn ~ 1 GeV is IP^at) 2 < 10 7 , cf. Table |n| Equation (33) then 


implies that Bi4 cff - ) < 10 12 for such decays. The proposed experiment at CERN-SPS (5DJ [HT] could produce D and 
D s mesons in numbers by several orders of magnitude higher than 10 12 , which would open the possibility to explore 
whether there is a production of sterile Majorana neutrinos N in the mass range Mn ~ 1 GeV. 

If the decays Bf —» /u. ± yu =1= 7r^ are consider ed (w e do not use P* decays as they are CKM-suppressed compared to 
P)f), the results of Figure [7jff and Equation (31b) imply an effective branching ratio 

(34) 


Br e ff(P^ —> 


10 2 |P 


fiN I 


which is similar to the case of D s , Equation (33) in a detector of the same length (L = 1 m, in our example). Since D s 


is significantly lighter than P c , the relevant neutrinos that give a sizeable effect are also lighter and thus longer living, 
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FIG. 4: The effective canonical branching ratio (32 I for the decays K ^ 
neutrino mass Mn- 






e, fi) as a function of the Majorana 




/Ww(GeV) /Ww(GeV) 


FIG. 5: The same as in Figure [ 4 ] but now for the decay of (a) D± mesons; (b) mesons. The solid lines are for £ = e, and 
the dashed lines for £ = (i. 


implying a smaller Pjy factor for the same detector length. On the other hand, D s decays have no CKM suppression 
compared to B c (|V CS | ss 1 while V c i, « 0.04) and the B c channels are more numerous, so that the true branching 
ratios of the latter are smaller. These two effects compensate in a given detector, as shown in Equations (33) and 
(34). However for a longer detector the observable D s branching ratio increases considerably mm- 


III. CHARGED PION DECAYS MEDIATED BY ON-SHELL MASSIVE NEUTRINOS 

In the previous Section we considered semileptonic decays of mesons heavier than the pion. For pion decays, there 
are purely leptonic modes only, since the pion is the lightest meson. In this Section we present and discuss the 
branching ratios for the LNC decay 7r =l= — > e ± N -A e ± e ± n T i' e and the LNV decay 7r =l= -a e ± IV -A e ± e ± We also 

consider the differential branching ratios dBr/dE M , where E^ is the energy of the produced /rC In contrast to the 
previous Section where the intermediate N neutrino has to be Majorana, here N can be either Majorana or Dirac. 
In the Dirac case, only the LNC mode is possible, while both LNC and LNV modes occur in the Majorana case. 
However, the experimental distinction between these two modes cannot be resolved by simply examining the final 
state, since the produced neutrino ( v e or v^) is not detectable. The major part of this Section refers to Ref. E3, and 
for certain details we use here results of Refs. |55J SS]. The formalism is somewhat more complicated now because 
we have four final particles (in the previous Section there were three). Nonetheless, several features turn out to be 
similar as in the previous Section. The considered processes are presented in Figures [8] and [9] As in the previous 
Section, we consider a scenario with at least one heavy sterile neutrino N, which has suppressed heavy-light mixing 
coefficients B^y with the first three neutrino flavors vg [l = e, /i, r), cf. Equation ([l]). The mentioned decay rates may 
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FIG. 6: (a) The effective canonical branching ratio ( |32[ ) as a function of Mjv for the lepton number violating (LNV) decays 
B ± —> y, w here £ = e,n (no visible difference between £ = e and l = [i ); (b) the theoretical canonical branching ratio 

Br, Equations (191 and (201, for these decays. 




FIG. 7: The same as in Figure [hj but for the LNV decays of the charmed mesons . 


be nonnegligible only if the intermediate neutrino N is on-shell, i.e., 

(M m + M e ) < M n < (M*. - Me) (35) 

and the process is of the s-type, Figures [H] and [9j Specifically, this means 106.2 MeV < Mjv < 139 MeV. 




FIG. 8: The lepton number conserving (LNC) process 7r + —¥ e + N —► e + fi e + v e which can be mediated by a Dirac or a 
Majorana on-shell neutrino N : (a) the direct ( D ) channel; (b) the crossed (C) channel. 
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FIG. 9: The lepton number violating (LNV) process 7r + —> e + N —> e + e + p which can be mediated by a neutrino N only if 
N is a Majorana particle: (a) the direct (D) channel; (b) the crossed (C) channel. 

A. Branching Ratios for —> e ± e ± p^v 

The decay widths T'( x ) (/T^ — > e ± e ± /i =F i/) (X=LNC, LNV) can be written in terms of the corresponding reduced 
decay amplitudes 


r (x) (7; 


e ± e ± p T v) = 


1 1 


2! 2(2tt) 8 


d4 \t^2\ 2 


(36) 


Here, 1/2! represents the symmetry factor from two final state electrons, and d 4 is the integration element of the 
phase space of the four final particles 


*= n 


\i =i 


d 3 Pj \ d 3 p/ d 3 p v 
2 E e {pj) ) 2E ll (p ll ) 21 p„ 


<5 (4) (Ptt ~Pi—P2-Pn~ Pv) 


(37) 


Here, p\ and P 2 are the momenta of e; in the direct channel, the e momentum at the first (left-hand) vertex is p\, and 
in the crossed channel it is P 2 , cf. Figures [i| and [ 9 J When we use the expressions for the amplitudes 7^± of Appendix 


(38) 


4, for the specific considered case N\ = A(and no N 2 ), the decay width ( |TTg| ) can be written as 

r^Tr* —► e ± e ± p^v) = |fe( x )| 2 [fW(i)I)*)+fW(C(7*)+f[ x 2(DC*)+f^ x ] ; (C'T>*] 

where X = LNC, LNV; are the corresponding mixing factors 

fclLNV) = B 2 n s /c(LNC) = B eN B* N 

and Hr (YZ*) are the normalized (i.e., without explicit mixing dependence) decay widths 


= Kl 


1 1 


1 


2! 2 AU (2?r) 8 


d 4 P (X) (F)P (X) (Z)* T^%(YZ*) 


(39) 


(40) 


The expressions for Tj2±(YZ*) are given in Appendix 4, for the direct (YZ* = DD*), crossed (Y Z* = CC*) 


and direct-crossed interference contributions (YZ* = DC*, CD*). We have T^)(DD*) = T) A /(flfl*) and 


AX), 


T^2(CC*) = tI X 1(CC*), hence the terms T^iDD*) and rl X ^(CC*) in Equation (381 have no subscripts ±. Further, 


in Equation ( |40| , p( x )(F) represent the N propagator functions of the direct and crossed channels (Y = D,C) 

1 


p(LNC) ^B) = 
p(LNC) (C) = 


[(Ptt - Pi) 2 ~ M% + iT N M N ] ’ 
1 

[{p«-P2) 2 -M% + ir N M N y 


p(LNV)( jD ) = M N p( LNG \D) 
P ( LNV)( C -) = M N P^ C) (C) 


(41a) 

(41b) 


where Tn = T(N —> all), and K 3 is the following constant: 




2.983 x 10~ 22 GeV 


-6 


(42) 
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It turns out that in the case when the intermediate N neutrino is on-shell, i.e ., when its mass is in the interval of 
Equation (35), the squares of the propagators (41) reduce to simple delta functions due to the inequality IV <C M,v 

2 


|P(LNC)(X)| 2 = 


1 


(Pn - Pk) 2 - Mh + iT N M N 


H(Pn — Pk) 2 — Ml) (IV < Mjy) 


AVIV 


(43) 


where pk = Pi,P 2 for X = D,C. In this on-shell case, the DD* and CC* contributions in Equation ( |38[ ) are large 
and equal, and the interference contributions DC* and CD* are negligible in comparison; we refer to |38| for details 
on this point. Hence the decay width (l38|) can be written in the on-shell case as 


e ± e ± pXv) = 2|fc (x) | 2 r(7T d 


0 ± P ±.,T 




(44) 


Here, the normalized decay width T^\DD*) = T^\CC*) = r(7T ± —> e ± e ± pXv) turns out to be the same for X = 
LNC and X = LNV, 


r(7T ± e ± e ± ^v) = r W{DD*) = T^(CC*) 

Kl Ml/ ,1/2/ 

= 192(27 t) 4 M 3 Tn^ \ x m Xe ^ ~ ^ Xe \ x K T 2 — x e )\ J- (x M , x e ) 
where the following notations are used: 

Kvu 2/2, 2/3) = 2/i + 2/1 + vl ~ 22/i2/2 - 2^2/3 - 21/32/1 


Ml 

Ml ’ 


M 2 




(45) 

(46a) 

(46b) 


and the function D(xij,,x e ) is given in Appendix 5. In the approximation M e = 0, the above result becomes simpler 


lim r(7r d 
M e -X) 


e ± e ± pXv) = 


Kl 


Ml 1 


192(2tt) 4 IVM3 


( x * - l) 2 /(v) 


where the function /( x^) — T(x^, 0) is 


/Vp) = 1 - 8.t m + 8x 3 - x A - 12a: 2 In 


(47) 

(48) 

(49) 


It can be checked that the decay rate (44), with N on shell, coincides with the factorized expression 

r(7T+ —> e + e + p,~v) = r(7T+ —> e + N)Br(N —$■ e + pTv) 

In Appendix 5 we also provide the differential decay rates dT^('!r ± —> e ± e ± Mv)/dE ll with respect to the final 
muon energy in the rest frame of N neutrino. They turn out to have quite different forms for X = LNC and X = 
LNV cases (we will return to this later in this Section). 

In order to obtain the branching fractions of the considered processes, we need to divide the decay width by the 
total decay width of the charged pion, r(7T ± —► all) 


r(7r + -> all) = 2.529 x 10" 17 GeV 


—G|/ 2 M 2 M ot \v u M 


M 2 \ 

1 ~m) {1+69n) 


(50a) 

(50b) 


where the expression (50b) represents the by far most dominant decay mode 7r d 

Ml (1 - Ml/Ml) 2 


pXi'n, and 


fig* = 


mi (i - m/M*y 


(51) 


represents a (very small) relative correction coming from the 7r ± — > e,u e decay ( Sg T ss 1.3 x 10 4 ). 

As we can see in Equations (44) and (451, the decay width r(7T ± —> e ± e ± p, T u) and its normalized counterpart 
TV 1 * 1 —> e^e^pXv) are inversely proportional to the (very small) decay width F(N —> all) = IV, which in turn is 
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proportional to the mixings 1C ~ \B^\ 2 (l = e, u, r), c/. Equations (16|-(|18[) and Figured] As in Section IIA 
effect represents the iV-on-shell effect Equation (43) and it makes the considered width Flyr* 1 


this 


e ± e ± // =F i/) by many 


orders of magnitude larger than it would be in the case of off-shell N. In the (narrow) mass interval (35) for on-shell 
N in the considered pion decays, the factor 1C in the width Tjy Equation (|16j) has the following approximate form 
(c/. Appendix 3 and Figure [2]): 




e ± e ± ^v) 


1.6\B eN \ 2 + 1.1(\B imN \ 2 + \B tN \ 2 ) 


(52) 


This expression, within the precision given here, is valid equally for the Majorana and for the Dirac N and agrees 
with that given in Ref. (38] for the Majorana case. However, the affirmation in Ref. [38] that 1C for Dirac N is smaller 
by a factor of two is not correct. 

In certain analogy with Section 
branching ratio Br^- 1 = r(7T ± —» 
dependence (and is independent of X = LNC or LNV) 


IIA 


we can define here the canonical branching ratio Biv as the part of the 
e ± /r =F i')/r(7r —> all) which contains no explicit or implicit heavy-light mixing 


Biv = 


2^rV 

-L 7T 


e ± e ± n T v) = 


1C 


|fc( x )|' 


( x )^± 


r Br^( 




li T u) 


1 K$M* 

16t T G|T(tt+ -> all) xl 
1 1 


2 xn(x 7r X/j,) 2 ( 1 4- 


A ^ (^7T5 Ip %e) [^7r 1 4~ ^e(*^7r 4“ 2 3?e)] *^"(*£/rp Xg') 

A ^ [Xm 1, X e ) [x^ 1 4“ X e (Xit -(-2 X e )] .A(x^, Xe) 


(53a) 

(53b) 

(53c) 


where the notations ( |16[ ) are used, and in Equation ( |53c[ ) we used the expression (125) in Appendix 5. We note that 
1C ~ \B tN \ 2 and |fc (x) | 2 ~ |FW| 4 , and BiA X) oc |A;( x > | 2 //C ~ \B eN \ 2 , where B lN stands for a generic heavy-light mixing 
coefficient {\B^ N \ ~ \B eN \ ~ |-B m jv|)- Naively, we should have very strong suppression Br^ a |fc( x )| 2 a \B( N \ 4 - 
however, the on-shellness (43) brings in factor 1 /Tjv oc 1//C ~ 1/\B(n\ 2 , reducing the suppression to Br^ oc |F^at| 2 . 
In Figure |l0| we present the canonical branching fraction Br^ as a function of M)v in the entire interval of on- 
shellness (|35b. In Figure 11 this canonical branching fraction is presented at the lower edge of the on-shellness 


interval, where the effects of M e ^ 0 turn out to be appreciable. We can see that the branching ratio is the largest 
when Mn ss 0.130 GeV. The differential branching ratios dBr^ /dE^, where E M is the final muon energy in the 



FIG. 10: The canonical branching ratio Br^-, Equation (531, as a function of the mass Mjy. The full formula was used (with 
M e = 0.511 x 10 -3 GeV). The formula for M e = 0 case gives a line which is in this Figure indistinguishable from the depicted 
line. 

V-rest frame, are obtained directly from the differential branching ratios dT^ /dE„, and the latter quantity for 
X = LNV is given explicitly in Appendix 5. The canonical differential branching ratios, free of any 
mixing dependence, can be defined in analogy with the definition of the canonical total branching ra¬ 
tio (53a), and, in contrast to canonical branching ratios (53) they do depend on whether X = LNC or 
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FIG. 11: The canonical branching ratio Br^ near the lower end point + M e (=0.1062 GeV): (a) in the interval below 0.107 
GeV; (b) in the interval below 0.110 GeV. The dashed line is for M e = 0, the full line includes the effects of M e = 0.511 x 10 -3 
GeV. 


X = LNV 


dBr^ — 2 c ^( 7r± e±e ~~ 

dE,, IV dE,, 


z H T v) 1C dBr (x ^(7 


e^e- 


z ^v) 


| fc( x )| 2 


dE„ 


(54) 


Explicit expressions for these quantities, when X = LNV and X = LNC, are given in Appendix 5 in Equations (126) 
and (128), and in the limit M e = 0 in Equation (129). 


The differential (and full) branching ratios for the process w ± —► e ± e ± differ in the cases when the on-shell N 

is Majorana or Dirac. When N is Dirac, only the X = LNC process contributes. When N is Majorana, both LNC 
and LNV processes contribute. We can write the differential and full branching ratios in these cases in terms of their 
canonical counterparts, by defining first the combined canonical differential branching ratios 


dBi'n^a) 

dE u 


/FW(LNV) 

dE u 


+ (1 - a) 


/FS - (LNC) 

rfBl V 

dE a 


(55) 


where 0 < a < 1. Then it is straightforward to check that in the cases of Dirac and Majorana N neutrino, the 
branching ratios can be expressed in terms of the above quantites (|55|) 


dBr ( - Dn '' ) (7r ± — > |^.(lnc)| 2 dBr w (a = 0) 

dE^ £( D ir.) dE^ 

dBr (Maj - ) (^ ± ->e ± e ± ^v) _ (|fc( LNV )|2 + |fc(LNC)| 2 ) rf Br OT (q M ) 

dE^ ^(Maj.) dEft 


(56a) 

(56b) 


where we recall the definition (39) of the coefficients and the “Majorana LNV admixture” parameter a« 

appearing in Equation (56b) is defined as 


| fc (LNV)|2 

~ (|jfc(LNV)|2 + |fc(LNC)|2) 


\B eN \ 2 

(\B eN \ 2 + I^nI 2 ) 


Integration of the relations (56) over leads to the full branching ratios 


Br (Dir ) (7T ± -» e ± e ± iX*v) 
Br (Maj ' ) (7T ± -T e ± e ± ^v) 


|fc (LNC) ! 2 p- _ \B eN \ 2 \B , N \ 2 p- 

K (Dir) " E^.^iTWl 2 ” 

(|fc( LNV )| 2 + |fc( LXC )| 2 ) \B eN \\\B eN \ 2 + \B, N \ 2 ) 


(57) 


(58a) 

(58b) 
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where we used the fact that the E^-integrated expression Br w (a) is independent of a , cf. Equations (53). In Figures 
we present these canonical differential branching ratios as a function of the muon energy E M in the rest frame of N, for 
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four different values of mass Mn in the on-shell interval (35), and for five different values of the admixture parameter 


(«m = 1, 0.8, 0.5, 0.2 and 0), The value olm = 0 corresponds to the case of Dirac N. From the curves of Figure 12 
we can conclude: if N is Major ana neutrino with a significant value of the admixture parameter olm and with mass 
in the interval (35), then the measurement of such differential branching ratios may be able to confirm the Majorana 


nature of N. If the differential branching ratios are studied with respect to the muon energy E in the pion rest 
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FIG. 12: The canonical differential b ranc hing ratio dBr n (oiM)/dE fl , as a function of the muon energy in the neutrino N rest 
frame, E^, as defined via Equations (55 I, for the decays mediated by a Majorana neutrino N, for various 


neutrino masses: (a) Mn = 0.112 GeV; (b) Mn = 0.119 GeV; (c) Mn = 0.126 GeV; (d) Mjv = 0.133 GeV. In each graph 
there are five curves, corresponding to different values of the admixture parameter om [Equation ( |57| )]: om = 1.0 is the solid 
(M) curve; 0.8 (dotted); 0.5 (dot-dashed); 0.2 (dashed). The case mediated by a Dirac neutrino (om = 0) is also presented as 
the solid line labelled (D). In comparison with Fig. 5 of Ref. [37] (where M e = 0 was taken and linear j/-scale was used), the 
true value of M e = 0.511 MeV is taken here and logarithmic j/-scale is used. 


frame, the distinction between the Dirac and Majorana case is more difficult, cf. Ref. E3- 


B. Effect of the Long Neutrino Lifetime on the Observability of ir* 


e ± e ± fPu 


The branching ratios presented in this Section so far, Equations (581 in conjunction with Equations (53), should be 
multiplied by the probability Pn for the on-shell neutrino N to decay within the detector (of length L), as explained 
in Section IIB Equations (22)—(26) and Figure[3] As a consequence, the effective (true, measurable) branching ratios 
are not Br (~ \Bim\ 2 ) of Equations (|58|), but Br e g = P/yBr 


Br 


eff 


) (tt ± -► e ± e ± n T iy) = P^ Dir - ) Br (Dir - ) ( 7 r 


N 


— ) /c (Dir - ) 

lm 


e ± e ± n T u) 
\B eN \ 2 \B, N \\- 


/c( Dir -) 


L \- - — 


= \B eN Y\B^ N \ 2 — P N Br n = \B eN \ 2 \B^ N \ 2 — Br Wi 


lm ) 




lm J 


•,eff 


(59a) 

(59b) 
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Br^ aj - ) (7r ± ->■ e ± e ± fjL T v) = P^ Maj °Br (Maj - ) ( tt : 


e ± e ± n T v) 


P, 


\ /c(Ma-j-) 


N \lm J 

= \B e N\ 2 (\B e N \ 2 + |-B/iJv| 2 ) 


|Peiv| 2 (|Peiv | 5 


IS, 


MiVl 


-Br* 


/C(Maj-) 

— J PjvBr*. = \B eN \ 2 (\B eN \~ 


I B 


| 


lm) Br7r ’ e * 


(60a) 

(60b) 


In Equ ations (59a) and (60a) we used the expressions (58) for Br and Equation (26) for Pjy In Equations (59b) 


and (60b) we introduced canonical (i.e., without any mixing dependence) effective branching ratio Br* je ff 

Br^eff = P jvBr* 

with Br* given in Equations (53), and the canoni cal nonsurvival probability Pn presented in Figure [3] in Section 


(61) 


IIB 


for a wide range of N neutrino masses. In Figure 13 we present Pjv in the here relevant narrower mass interval (351. 



A4(GeV) 


FIG. 13: The canonical proba bility Pn = (lm)/Ljv, Equation (261, as a function of the neutrino mass A-f.y in the approximate 
interval where it is on shell (35 1 for decays n ± 

[= (l-/3y)“ 1/2 ] =2. 


e ± e ± /r =F //. The Lorentz lab time dilation factor is chosen to be 7 jv 


As in the case of semileptonic LNV decays of Section II B[ we notice that also here in the effective branching 
ratios, Equations (59) and (60), the com plicated dependence on the mixing parameters entailed by the factors K, = 
^2gAT in I B(n 1 2 of iv [cf. Equation (18)], cancels out, and there remains only simple dependence on the mixing 
parameters, in the form \B s n\ 2 \B^n\ 2 or |P e Ar| 2 (|P e jv| 2 + I-B^atI 2 ). Further, in comparison with the branching ratios 


Br of Equations (|58|), which are ~ |P^at| 2 , the effective (true) branching ratios Br e ff of Equations (59) and (601 are 


unfortunately significantly more suppressed by the mixing parameters, namely Br e ff ~ |P^at| 4 - 

The prese ntly known experimental bounds on the mixing parameters \B^\ 2 (£ = e, /it, r) in the here relevant narrow 
mass range (35), are: |P e Af| 2 1$ 10 -8 [§5]; I-B^tvI 2 1$ 10 -6 (100H102) : |P r Ar| 2 < 10 -4 |103| : cf. also Refs. [!M 1711 [IM1 


The future pion factories, such as the Project X at Fermilab, will be designed to produce charged pions with lab 
energies E * of a few GeV and luminosities ~ 10 22 cm _2 s _1 [1(7B , 1107 !. and ~ 10 29 charged pions could be expected 
per year. 


The canonical effective branching ratio (61) can be estimated as 


Br 


7r,eff 


< 10 


-6 


(62) 


as can be inferred from Equation (61) and Figures [13] and ITOl Equations (59) and (60) then imply that the effective 


(true) branching ratios for the considered reactions are (we assume L = 1 m for the detector length) 


Brig 1 ' 0 < \B eN \ 2 \B^ N \ 2 10~ 6 , 

l^eiv| 2 (|Reiv| 2 + |B^| 2 )10- 6 


(63a) 

(63b) 
















































17 


If the larger among the mixing elements (|£^at| 2 , i = e,n) is l-B^jvl 2 (< 10 6 ), the LNC processes dominate, and 
the effective branching ratios (63) have the common upper bounds 


Br 


(Dir.,Maj.) 

eff 


< I b, 


eN\ 2 \B 


)iN 


| 2 1CT 6 < IP^PHT 12 


(64) 

This 


If in this case |P e jv| 2 is close to its present upper bound, \B £ n\ 2 ~ 10 —8 , we obtain Br^' r ' : aj '' < 10 _2 ° 
implies that up to 10 9 events 7r ± —> e^e^yffv could be detected per year in such a scenario. 

On the other hand, if the larger among the mixing elements (B(n\ 2 , £ = e, ff) is |P e iv| 2 (< 10 —8 ), the LNV processes 
dominate, and the effective branching ratios (63) have the following upper bounds: 


Br 

Br 


(Dir.) 

eff 

(Maj.) 

eff 


< 


S„wl 2 10- 14 


l*N\ 
< 10" 22 


(65a) 

(65b) 


In such a case we have |P m at| 2 < |P e iv| 2 < 10 8 , and up to 10 7 events could be detected per year. 
The present upper bounds on \B^n\ 2 suggest that the first scenario, Equation (641, is more likely. 


The measurement of the effective branching ratios alone cannot distinguish between the Dirac and the Majorana 
character of intermediate neutrino N. However, as argued in Section |III A| and presented in Figure [l2j the mea¬ 
surement of the differential branching ratios for the considered processes is a promising way to discern the character 
of the neutrinos. However, the differential branching ratios (56) must be multiplied by the nonsurvival probability 


P/v, in order to obtain the effective (true, measurable) differential branching ratios dBr^/cLE^. In analogy with 
Equations (591 and (60) we obtain, using Equation (56) 


dB r£ ir 'V 


e ± e ± 




dE„ 


= Pi 


dBr( Dir -)(7T 4 


e ± e ± 


N- 




= I B, 


eJV 


dEp 

2 11? | 2 ( B ^ dBr 7rje ff(o: = 0) 


lm/ 




,=t„± 


!X*v) 


dE u 


= P 


N~ 


dBr (Maj - } (7r d 


0 ± 0 ± 


dE^ 


dEu 


\B eN \ 2 (\B eN \ 2 


T> |2\ { \ dBY-K 5 eff(o:M) 

B '“ V| > (to) dE, 


(66a) 


(66b) 


where the canonical differential effective branching ratios were introduced, in analogy with Equation (55) 


^Br 7r)eff (a) _ ^ dBr n (a) 

-,,,-= — e -— 


dE u 


dEu 


(67) 


12 


with dBv n (a) / dE^ given in Equation (55). 
tion (57). In an alog y with Figure 
y -axes of Figure 


The parameter a« appearing in Equation (66b) was defined in Equa- 


and using the values of Pjv from Figure 
12h-d, must be multiplied by Pat ~ 2.0 x 10” 3 , 2.8 x 10 -3 
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we deduce that the values of the 
3.7 x 10” 3 , and 4.8 x 10~ 3 , respectively 


(assuming L = 1 m and 7 jv = 2), 
effective ratios dBr W)e a(a)/djB M . 


in order to obtain the representation for the curves of the canonical differential 


IV. CP VIOLATION IN CHARGED MESON DECAYS MEDIATED BY MASSIVE STERILE 

NEUTRINOS 

CP violation in the lepton sector could be measured by neutrino oscillations I66| . Here we consider the possibilities 
of measuring CP violation in meson decays mediated by sterile neutrinos N , such as the semileptonic LNV decays of 
charged heavy pseudoscalar mesons considered in Section [ilj or the leptonic (LNV and LNC) decays of charged pions 
considered in Section Hill 

It turns out that CP violation in all such decays is possible in scenarios with at least two massive sterile neutrinos 
Nj ( j = 1,2). CP violation in the neutrino sector is expected whether neutrinos are Dirac or Majorana particles. 
However, in the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix jl08Hll0j . the number of possible CP- 
violating phases is larger when the neutrinos are Majorana particles. If n = 3 + Af is the total number of neutrinos 
(A f is the number of sterile neutrinos), the number of CP-violating phases is n{n — l)/2 if neutrinos are Majorana, 
and (n — l)(n — 2)/2 if neutrinos are Dirac, cf. Ref. Ill'l l. 
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CP violation in the decays 7r =t —► e ± Nj — ► e ± e ± p/ f v was investigated in Ref. j35], and in the decays M ± -a £- t Nj -a 
^ in Refs. Esi mi- In both cases, it turns out that, even though the rates are extremely small, the CP violation 
asymmetry in these charged decays may become appreciable, even close to order unity, when the two intermediate 
neutrinos can go on shell and are almost degenerate in mass Mm, « This is to be contrasted with CP violation 

in charged meson decays due to the standard CKM mechanism in the quarks sector, where the rates are larger but 
the asymmetries are much smaller, e.g., of order 10~ 4 in K ± -A 37 t decays |112j . 


A. CP Violation in Semileptonic LNV Decays —> iflfM ~ 

As mentioned above, we will consider the scenario with at least two sterile neutrinos, Nj (j = 1,2), both of which 
can go on shell in the intermediate state, i. e.. with masses satisfying the condition shown in Equation (l2|). The 
processes of interest are again those of Figure [l] except that now in both the direct (D) and crossed (C) channels 
there are two possible neutrinos exchanged: Ni or .ZV 2 . The relative measure of CP violation for these processes will 
be the asymmetry: 


Acp(M) = r(M 


- T(M+ -A £f£+M'~) 


T(M~ 


^qM'+) + r(M+ —> i+i+M'-) 


( 68 ) 


The corresponding LNV decay widths T(M ± -a T ) will be obtained now, in the scenario of two sterile 

neutrinos (A f = 2), in close analogy with the calculation in Section IIA which was performed for the case of one 
neutrino N (A f = 1). The relations |3| and © and |9]) there hold without change now, but the relations © © 
obtain the following slightly more general form when AT = 2: 


T(M ± -A ixtiM T ) = 

2 2 

[r{DD*)ij + T(CC*)ij + P ±{DC*)ij + T±(CD*)ij] 


(69) 


*=1 j'=i 


Here, are the mixing coefficients 


k { ~ ] = B (lNj Be 2Nj , k { + ] = (kj->)* (70) 

and r±(XV*)ij are 2x2 matrices, and represent the normalized (i.e., without the explicit mixing dependence) 
contributions of Nj exchange in the X channel and complex-conjugate of the Nj exchange in the Y channel (X, Y = 
C, D) 


r±{XY*)ij = K 1 


1 1 


1 


2! 2 M m (2tt) 5 


d 3 Pi{X)Pj ( Y)*M Ni M n T±(X)T±(Y)* 


(71) 


The expressions for T±(X)T±(Y)* (where X, Y = D,C) are the same as in Equation Q and are given in Appendix 
1, and Pj(X) (X = D , C) are the propagators of the exchanged neutrinos Nj in the direct and crossed channels 


PAD) 


PAG) 


1 

(pm ~ Pi) 2 ~ + iT Nj M Nj 

1 

( Pm - P 2) 2 ~ + iT N .M N . 


(72a) 

(72b) 


We will disregard effects due to non-diagonal neutrino widths in their mass basis. For these details we refer to 


Ref. 


The total decay width of Nj, T N 


parameter AA Equation 


1C. 


— A/jyv. 


is given by Equations (16) (18), where each Nj has its own mixing 
'Bj N ,| 2 , where the coefficients A fw- as a function of M N . are given 


by the left-hand Figure [2] in Section II A| 

As we will see below, the CP asymmetry parameter Acp(M), Equation (68), 
nificant nonzero value if simultaneously: (a) the phases <f>ij of the PMNS heavy-light 


may acquire a sig- 
mixing elements 


B(.Nj = \BiN j \ exp(i(f>ij) fulfill certain conditions: |sin0 2 i| = | sin(</>^ l2 + <j>t 2 2 — (fit, 1 — <f>e 2 i)\ A' 1 ; and (b) the 
mass difference A Mn = Mn 2 — Mn, is sufficiently small (|AMy| . 
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First let us calculate the quantities 

S T (M) = F(M~ -> =f r (M+ -»■ £+£+M'~) (73) 

appearing in the numerator and the denominator of the CP violation parameter Acp(M) Equation ( |68| . We introduce 
the following notations which will be needed below: 


_ \Bt lN2 \ _ \B t2N ,\ 

\BcinA ’ t2 I B t2Nl \ 

Be k N j = | B(, kN .\e l<t>t ^ {k,j = 1,2) 

@ij = (&ii d“ 0^2& 0 b,7 0^2j) (bj 1; 2) 


(74a) 

(74b) 

(74c) 


For e xamp le, in the specific case when £i = £2 = £b we have $21 = 2(^2 — 0/xi) = 2(arg(£? Al jv 2 ) — arg(S /i Ar 1 )). As in 
Section IIA when both A(j are on-shell, it turns out that the interference contributions to the quantities Szp(M) from 
the direct ( D ) and crossed (C) channels (DC* and CD*) are suppressed by several orders of magnitude in comparison 
with the contributions from the direct (DD*) and crossed (CC*) channels, and we will neglect them (we refer to [35] 


for details on this point). Then it follows from the expression (69) 


S-(M) = F(M~ £^£~m'+) - F(M+ -)• £+£+M“) 

= 4(2 - Se 1 e 2 )\Be lN X\Be 2Nl \\Be 1 N 2 \\Be 2 N 2 \{sin9 21 [imf (DD*) 12 + ImF^C*)^ } 


(75) 


and 


S+(M) = F(M~ -»• M'+) +T(M+ ->■ £+£+m'~) 

121 


= 2(2 - S ili2 )\B eiNl \ 2 \Be 2Nl \ 2 < r (DD*) U 




1 + Kg. 


2 2 T(DD*) 22 


T(DD* 


+ 2k^ 1 cos 02 i^i 


+r(C'C'*) n 


r(cc*) 2 2 

1 + ftf, ~ -b 2 K(, 1 K£ 0 COS 0 2 lSi 

1 2 r(CC'*)u 


In the sum (76), the coefficient <5i represents the effect of N- t -N 2 overlap contributions 

Ref (A'X *) 12 


Sj = 


T(XX*) 


(■X = D-C ; j = 1; 2) 


(76) 


(77) 




We expect <5i « 0 when A Mm F^ (where: A Mat = Mjv 2 — My, > 0); numerical calculations (see later) confirm 
this expectation and show that Sj is practically independent of the channel X = D , C. The normalized decay widths 
T(DD*)jj and F(CC*)jj are those of Equations (12) and (13) of Section IIA with the substitutions My K> M,y ; , 
Un eb y Nj = M 2 n ./M 2 m [cf. Equation (14b)] and F N i-b F Nj (= K.jF N ) 

f(A>£>*)„■ = ^-^f^ L X 1/2 (l,y Nj ,y£ 1 )X 1/2 ( 1 , —, —) Q{y Nj \y^yt 2 ,y') (78) 

128M V y Nj y Nj J 

where j = 1 or j = 2; r(C'C'*) J j is obtained from F(DD*)jj by the simple exchange y ^ o y i2 [cf. Equation (13)]. 

For evaluation of the CP-violating difference S-(M), Equation (751, the quantity Imr(XX*)i 2 
(X = D: C) is of central importance. In the integrand of Imr(XX*)i 2 appears as a factor the following combi¬ 
nation of the propagators of N\ and N 2 [cf. Equation ©]: 

( Vn ~ M % ri ) Fn 2 Mn 2 - Tv, M Nl ( p 2 N - M ^) 


Im (P 1 (D)P 2 (£>)*) = 


(p 2 n - m 2 Ni ) + f 2 n m 2 Ni \ [(p 2 n m 2 N2 ) + r 2 n m 2 N2 
V x ^2 __ M 2 — M% 2 ) + S(p 2 N - MyJ] 


In Equation (79b) we used the narrow width approximation: 

F Nj M Nj 


(p 2 n~M 2 Nj ) +F 2 n M 2 Nj 


— irS(p N — M n ) (for ryc AfjVj) 


(79a) 

(79b) 

(80) 
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and in Equation (79b I the parameter 77 was introduced which parametrizes any deviation from the naive 
expectation 77 = 1. We expect rj ~ 1 when AMj P/v, , Ta^ , where AAfjy = M^ 2 — > 

0. In Appendix 6 we argue that this parameter 77 , in the case of near degeneracy A M N <C M Nl , 
is a simple function of only one variable y = AM/Tn , where AMn = Mjv 2 — > 0 and 

Ttv = (1/2) (r^ + IV,,) 


v(y) = 


y 


( y 2 +1) 


tj=AMn/Fn 


when AM n <C M Ni = M N (=> AM^ = 2M N AM N ), and where 

_ 1 /-p , t- A - AM N 

Tat = -(rjVi + rjv 2 ), y = 


(81) 


(82) 


This implies that in the case of two almost degenerate sterile neutrinos (AAV "C AV,) we have for the factor in 
Equation (79b) the following identities: 


77 x 


1 

A M% 


1 v(y) 

2M n T n y 


AM 2 n 


(AM 2 N y + 4M 2 n T 2 n 


(83) 


We note that the factor 4 in the denominator on the right-hand side of Equation (83) is nontrivial, because a 


somewhat different result would have been obtained by a more simple and direct consideration of the expression (79a) 
in the limit IV M N . (j = 1,2). The result (81) [or equivalently, Equation 

numerical evaluation of ImT(A'X*)i 2 , Refs. 


has been confirmed also by 
(see below). The mechanism |79|) [with the identity |83| ] is of 
central importance for the CP violation in the processes considered here. The quantity 77 /AAf^-, at fixed Tat and 
fixed Mat = AV 1 , achieves its maximum when AM^ = 2MnTn, i.e., y = 1 [<t=> AMn =^n (<C AV)], Ae., when 
the two sterile neutrinos are almost degenerate. If AMn /IV (*-e., y / 1), then the quantity y(y)/y = y/(y 2 + 1) 
in Equation (831 is very small and CP violation effects disappear. 

in the context of CP violation of leptonic decays of charged pions, and 


The mechanism (79) was used in Ref. 


in Refs. HE SHI in the here presented context of CP violation of semileptonic LNV decays of heavy pseudoscalars. 


We recall that the expression (79) has the same structure with Dirac delta functions as Equation (10) in Section IIA 


however, the factors in front of these Dirac delta functions are different now. Therefore, integration over the final 
particles’ phase space can be performed now in the same way as in Section m ^ analytically. This leads, in 
analogy with Equation (12), to the result 


1 K 2 M 5 m 


Imr(DV)i 2 =„(„) x ^ 


X EA 1 / 2 ( 1 ,V.^ 1 )A 1/2 (l, —, — ) Q{yN j ;ye 1 ,ye 2 iy l ) 
fri V ay, ax, 1 


ImT(CC*) 12 = lmT(DD*) 12 ( yil ^ y i2 ) 


(84a) 

(84b) 


where we recall the notations Equations (14) and (82), yN = Af/./A7/ f . the function Q is presented in Appendix 2, 


and we denoted AMn = Mn 2 — Mn 1 > 0 and AMfr = M^ 2 — ATj^ > 0. We note that in Equation (84) we have not 
yet assumed the near degeneracy of the two sterile neutrinos. 

From here on in this Section, we will consider the case of near degeneracy of the two on-shell sterile neutrinos 
(AAV -C AV, where AV = AVi), in which case Equations (81) and (82) hold and the quantities Equations (83) 


and (|84|) become appreciable and CP violation can thus become significant. Therefore, we have 


vn 2 ~ y A?! = vn — 


Ml 

Mm 


(85) 


where M N = M Nl ~ AV 2 • In this case the identity (81) holds, cf. Appendix 6 , and the expression (84a) becomes 
simpler 


y{y) K 2 M 5 m M n 


lmT(DD*)i 2 = x 

v(y)^^r^^ 2/Cl 


y 1287 r 2 rAT 
/ . 

-r (DD m 


y 


X 1/2 (l,y N ,ye 1 ) x 1/2 fl,—,—) Q{yN',ye 1 ,ye 2 ,y') ( 86 a) 

V VN Vn J 

( 86 b) 


(/Cr+/C 2 ) 
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TABLE IV: Values of theAfi-A^ overlap parameter 5(y) as a function of y = AMjv/Tjv. 


_ A M n 

y- r„ N 

logio y 

<%) 

0.10 

-1.000 

0.989 ±0.001 

0.30 

-0.523 

0.917 ±0.001 

0.50 

-0.301 

0.800 ±0.001 

0.70 

-0.155 

0.673 ±0.001 

0.80 

-0.097 

0.610 ±0.001 

0.90 

-0.046 

0.551 ±0.002 

1.00 

0.000 

0.499 ± 0.002 

1.25 

0.097 

0.390 ±0.003 

1.67 

0.222 

0.264 ±0.003 

2.50 

0.398 

0.138 ±0.001 

5.00 

0.699 

0.038 ±0.001 

10.0 

1.000 

0.0098 ±0.0010 


where in the last identity we used the expression ( 12 ), the identity (82), and the fact that 


r jv 2 /rjVj = IC 2 /IC 1 [cf. Equation (18) for Ni and N 2 


The normalized decay matrix elements r±(AT*)y, Equation (71), were evaluated in Ref. [35] also numerically, by 


Monte Carlo integration and using finite small widths Ey in the propagators. The numerical calculations confirmed 
the presented formulas, among them the expressions (78), ([ 86 ]), and (f8T|)— (83). The form (81) of 77 ( 7 /) was confirmed 


numerically with a precision better than a few per mille. The numerical evaluations also confirmed that the direct- 
crossed interference terms ( DC * and CD*) are really negligible. We refer to [35] for details. 

Further, the mentioned numerical evaluations gave us values of the N\-N 2 overlap parameter <5i as defined in 
Equation ©, i-e., the parameter which appears in the expression ( f7(i| ) and represents the N^-N 2 overlap effects. 
It turned out that the numerical values of the parameters Sj (j = 1,2), as well as of 77 , are practically independent 
of: the channel contribution considered ( DD * or CC*), of the type of pseudoscalar mesons (M ± , M T ), and of 
the light leptons (£ 1,^2 = e,/r) involved in the considered decays. The numerical results show that the parameter 
S = (l/ 2 )(Ji ± <5 2 ) is a function of only one variable, namely y = A M^/T yr (the same is true for 77 ) 


<5 = 


< 5 2 


<5(7/), 5=^(S 1+ S 2 ) 

(87a) 

f (DD*) 22 _ r Nl _ fc t 
f{DD*) 11 Tjv, jc 2 

(87b) 


The numerical values of the parameter <5 are given in Table IV as a function of y. It is not clear whether there 
exists a simple analytic expression for <5 as a function of y. Further, in Figure 14 we present the quantities 77 /y and <5 
as a function of y. The values for <5 in Table |IV] are practically equal to the values of the corres ponding 6 parameter 
in the rare leptonic decays of the charged pions 7 i - 


--N ■ 


e ± e ± /i =F i 7 , cf. next Section |IV A| and Ref. 


Branching r atio s of experimental significance here can be defined by dividing the expressions S^(M), Equations (73) 
and (75) and (76), by the corresponding sum of total decay widths [T(M + —> all)+r(M~ —> all)], which is practically 
equal to S> all) 


Br+(M) 

Br_(M) 


S+(M) 


2 T(M± -> all) 
M CP (M)Br + (M) = 


S_(M) 
2T(M ± -)• all) 


( 88 a) 

( 88 b) 


If employing the c anon ical (independent of mixing) branching ratio Br(?/jv; ye 1 , ye 2 , y') = Bt(DD*) defined in Equa¬ 


tion (20) of Section 


IIA 


and its CC* analog, Br(C'C*) = Br (yN',yi 2 ,ye 1 ,y l ), then it turns out that the branching 
ratios Br±(M) of Equation ( 88 ) can be rewritten in terms of Br, of the heavy-light mixing parameters |£W | and 
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y 


y 


FIG. 14: The suppression factors r](y)/y and 5(y) as a function of y = AMjv/Tjv: (a) for y < 1 (on the linear y scale); (b) for 
1 < y < 10 (on the logarithmic y scale). 


K.j = Yh-^tNj\B(Nj | 2 , of function rj(y)/y = y/(y 2 + 1) and the overlap function S(y) tabulated in Table 


IV 


Br +(M) = 


t{m~ -+ m'+) + r(M+ -► e+i+M'-) 


2r(M± -H. all) 


= 2(2 -S tlt2 ) 


E 

1=1 


1-SfiJV,- \ 2 \Bt 2 Nj 


AC, 


+45 (y) \ B ^ n i\\B+nA\B+ N2 \\B 12N2I cqs02i 


(ACi+ACa) 


(Br(D£>*) + Br((7<7*)) 


Br_(M) = A C p(M)Br + (M) = 

= 8(2-^ 2 ) 


r(M- -> - r(M+ -> m'~) 


2T(M± -> all) 

^ 1 iV 1 ||^ 2 iV 1 ||^ 1 iv 2 ||B, 2 ^ 2 | t?(2/) 


(ACi+ACa) 


sin09i (Br(£»D*) + Br(<7<7*)) 

y 


(89a) 


(89b) 


This leads to an expression for the CP violation parameter ylcp (M) defined in Equation (68), which now involves 
only the heavy-light mixing parameters \Bin.\ and AT, [cf. Equation (18)], the function r/{y)/y = y/(y 2 + 1), where 


y = AMjv/T jv, and the overlap function 6(y) tabulated in Table IV 


A (M) _ Br_(M) T(M~ -> 1^2 ^' + ) - r(M+ -> l+£+M'~) 

CPl ’ Br+(M) T(M~ -> + T(M+ -> (.+(.+M'~) 

4 sin 02i y 


Ei-irg-- (/Cig-zc^) _|_ 4j(y) cos0 21 j (y 2 + 1 ) 

1 \b 1iNi \\b 12Ni \\b 1iN2 ||rj£ 2 jv 2 1 

4 sin 0 2 i y 

{ K Ke *( 1 + ’k) + ^ ( x +§ 7 ) ] + 4< ^) cos 021 } (2/2 +1} 


(90a) 

(90b) 


In the usually considered case t\ = i 2 (= £), i.e., when the considered decays are M+ -+ £+£+M T , 
the formulas (89) and |90| get simpler, because in such a case 61(17(7*) = Br(DD*) = Br, and 
Be 2 Nj = Bi lNj = B( Nj , and K( 2 = kh x = Ke 


Br +(M) = 4 


+ 4 <%)cos 0 21 


i=l ^ 


(/Ci+/C 2 ) 


Br_(M) = ^ C p(M)Br+(M) = 16 ^S!-L^E_ sin 0 21 T4f^ B r 


Br 


Acp(M) = 


(ACr+ACa) 

4sin0 2 i 


2/ 


2/ 


{b 2 (‘+fe)+^ (‘+ft)] + «(»>■=» fe} (b2 + 1} 


(91a) 

(91b) 

(91c) 
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These formulas become even simpler if the absolute values of the heavy-light mixings of Ni and N 2 are equal (but 
not their phases), i.e., when 


I B 


'e'No 


\ B e ' Nl \ 2 = \B 


e’N I 


= e,//, r) 


(92) 


In such a case, we have all « 1, and K .2 « /Ci = 1C, and therefore the expressions (91) reduce to 


Br+(M) « 8 - 


l£N\ 


-Br (1 + 0(5)) 


K 

Br_(M) = A C p(M)Bt+(M) w 8 ^|^- sin 6 » 2 i ^B? 
A C p(M) » sin0 21 -^(l + 0(5))< £sin 6 l 21 (1 + 0(5)) 

V + 1 2 


(93a) 

(93b) 

(93c) 


In these expressions, we assumed, in addition, that the Ni~N 2 overlap terms are small (0(5)). 

In order to obtain the corresponding effective (true) branching ratios, we have to multiply the above branching 
ratios Br± by the decay-within-the-detector probability P\r. = P N ICj(L/ lm), as in Section IIB (L is the length 


of the detector). Again, the complicated mixing dependence entailed in the parameters ICj gets cancelled in this 
multiplication. When adopting the simplifying assumption Equation ( [92] ) , i.e., the validity of Equations (93), we 
obtain the following effective branching ratio Br e g and the CP violation effective branching ratio AcpBr e fj: 


Br eff (M ± —> e ± e ± M T ) = P n Bt+(M) 


q 

1^ 

K. 

[ 8 |b “ |4 b;1 

V y 


K. 


= ( 8 |-EW| 4 PatBi- = 2\B ( N \ 4 ^ —^ Br e g 


AcpBr eff (M ± -+ e ± e ± M'*) = sin 0 2 \—^—— x 2\B m \ 4 ( AV ) B r eff 

y z + 1 \lm / 

,4 ( L 


< I B. 


'IN I 


l lm 


|4 ( L \ 

sin 0 2 iBr efr 


(94a) 


(94b) 


In these formulas, we used the canonical effective branching ratio Br e ff as defined via Equations (32) and (20) 
and depicted in Figures 00 as a function of Mjv (= M/v, ps Mjv 2 ). The values of the Lorentz factors in the lab 
system are taken to be = 2 for both Ni and iV 2 , keeping in mind that Br e g scales as I/ 7 JV. We recall that 
621 = 2 ((j)i 2 — fa 1 ) = 2(arg (Bin 2 ) — aig(BtN 1 )). We notice that on the right-hand side of Equation (|94a[) there is an 


additional factor two in comparison with Equation (31b) of Section IIB this factor two comes from the fact that we 
now have contributions of two intermediate neutrinos Aq and N 2 , and we neglected the contributions from the Ni-N 2 
overlap (0(5)). 

As at the end of Section 


IIB 


let us consider now as an illustrative example the decays Df 


let us assume that \B^jsr \ 2 is the dominant mixing. In such a case, the estimate Equation (33) is still valid, and the 


7T 1 


In addition, 


CP-violating difference of the effective branching ratios, AcpBr e ff, is obtained by comparison of Equations (94a) and 
(|94b| 


AcpBr ef i(T>s /P/i*7 t t ) ~ 10 2 |P M 7 v| sin6» 2 i 


y 


y 2 +1 


ij=AMn/Fn 


< 10 2 |P M Ar| 4 sin 0 2 \ 


(95) 


Since for Mn ps 1 GeV we have at present |P^at| 2 5, 10 


| 2 < 10~ 7 , cf. Table |n| Equation |95l ) means that AcpBr (cff ^ < 10 -12 
for such decays. As already mentioned at the end of Section IIB the proposed CERN-SPS experiment p3Qj[8l] could 
produce D and D s mesons in numbers by several orders of magnitude higher than 10 12 , and production of the 
sterile Majorana neutrinos Nj could be explored. Further, if there exist two almost degenerate sterile neutrinos of 
mass Mn ~ 1 GeV (this is so in the pMSM model (BO, Ell 1B7K72] !. such that y = AM,y /Yn ~ 1, then we would 
have r](y)/y = y/(y 2 + 1) ~ 1. In such a case the estimate ( |95| would imply that the CP-violating difference of 
effective branching ratios, AcpBr( e g) (D s ), is of the same order as the effective branching ratio Bi-( e «)(D a ) (if the 
phase difference |0 2 i| 1). Therefore, if experiments can discover the mentioned pMSM-type Majorana neutrinos, 

they will possibly detect also CP violation effects coming from the Majorana neutrinos. 
































24 


B. CP Violation in Pion Decays 7p 


e ± e ± // =F p 


In this Section, we will only briefly outline the calculation of t he CP violation asymmetry in the (LNC and LNV) 
semileptonic decays 7 r =l= — t e ± e ± n T i' as described in Section III We will assume the presence of at least two n early 
degenerate sterile neutrinos Nj (j = 1,2) that can go on shell in the intermediate stat e, a s in Section IV A The 
present Section is a similar extension of the analysis of the decays 7 r =l= — t e ± e ± ^ of Section III to two sterile neutrinos. 


The results of the present Section are largely based on Ref. [38]. Only few details will be presented here, for further 
details we refer to Ref. [35] ■ 

Similarly to the previous Section |IV A] the quantities relevant for the CP violation will be 


Br 


(X) 

7T,± 


S±P(n) _ P( x )(7r —t e e [i + v) ± r( x )(7r + —t e + e + /z v) 


A (x) 

CP 


2r(7T± 

Brgl 

Brg 


all) 2r(7r± -t all) 

r( x )(7T^ —t e~e~ /i + p) — r( x )(7T+ 

~ r( x )(7T- —t e~e~fi + iy) + r( x )(7T+ —t e+e+/z _ p) 


e + e + /i v) 


(96a) 

(96b) 


where X = LNC, LNV. The total branching ratios are Br± = Br^ NV ^ + Br£ jN< A when Nj are Majorana neutrinos, 

and Br± = Br^ NC ^ when Nj are Dirac neutrinos. We adopt the same conventions and the same notations as in the 
previous Section [IV A| In addition, since we have now LNV and LNC processes, we introduce the additional notations 


0 (LNV) = 2(</> e2 — < j > e i ) 

0 (LNC) = ^ e2 _ 0 el ) _ 2 - (j)^) 


(97a) 

(97b) 


As in Section IV A the requirement that the quantities sin 9^ (here: X = LNV, LNC) be nonzero, and the 
requirement of the near degeneracy of the two neutrinos (AMjv <C Mjv, = Mjv) in conjunction with the expressions 
Equations (79)—(83) for Im(Pi(P)P 2 (P)*), are needed in order that the CP violation parameters A\ c P 7 ^ 0 acquire 
nonnegligiblevalues. Analysis similar to that of the previous Section IV A (but algebraically more complicated) leads 


then to the results for the quantities defined in Equations (96). More specifically, the results for the Dirac case, 
Br^f -* = Br^ LXC) and are the following: 


Br 


(Dir.) _ P( LNC >(^ -t e~e~fi + v) + r( LNC )( 7 r+ e + e + n~v) 


7T,+ 


a: 


2r(j r± -t all) 


2 I2|p„,2 


E 

3=1 


\BeNj 




/C, 


' \BeN 1 \\B e N 2 \\B^ Nl \\B ^ N2 \ 0 (LNC) 

(/Ci + /C 2 ) 


Br, 


(98a) 


(Dir.) _ r (LNC) (7T —tee /i + i/) — r( LNC )(7r+ —t e + e + fi v) 
7T, CP 


rtLNC)^- 


e~e~fi + iy) + r( LNC )(7r+ —t e + e + ^~v) 
4 sin 0( LNC ) 


v(y) 


\B cNl 1 | B^ Nl 1 

(l , K.2 \ 

_|BeJV 2 1 \B 11 N 2 1 

V 1 + xJ 


pw pw ( x + &) + 4S ^ cos0(LNC) 


|rv2 1 


(98b) 


The expression for the canonical quantity Br„, appearing in Equation (98a), is given in Equa¬ 


tion (53) in conjunction with the notation (46) in Section III A The results for the Majorana case, 
B r (Maj. b 1 C lxv 1 _|— Br ^. LXC ' 1 and A^?p^ are the following: 


Br 


(Maj.) 

7T,+ 


E x = 


X=LNV,LNC 


(r^ x ^(7r —tee /z + i/) + r( x )(7r + 


e + e + ji z/)) 


2r(7r 


r 2 

E 

3 =1 


B, 


eNi 


I B 


eN t 


IP 


fiNj 


/c, 


all) 

-4 S(y) 


IP, 


eN\ | |PeAi 2 


(AC 1 +AC 2 ) 


x (\B eNl \\B eN2 \ cos 6» (lnv) + 


B 11 N\ 11 P/r N 2 I COS 


0(LNC)\ 


Br, 


(99a) 
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^ (Maj.) 

•^TT, CP 


X^X— lnv.lnc ( r(x) ( 7r e ^+i/)-rW(7T+ 

Sx=LNV,LNC (r^ X) ( 7r —> e e H + v) + r fX ^(7T+ 

4 fsin^ LNV ) + l^ill^l sin0 (LNC)X 
V \B eNl \\B eN2 \ ) 


—> e + e + fi is)) 
—> e + e + y~is)) 


C|-B e JVi | 


I B 


M ATi| 


\BeN 2 \ 2 


+4(5(y) (cos0 (LNV) 


IS 


i + & 

/iNi 11 S^JV2 


(|S e AT 2 | 2 + |S m at 2 | 2 ) 


|S eJVl | 2 

-1 


\BeN\ | |S e jv 2 


COS 


0(LNC)j 


v (y) 



(99b) 


The function y{y)/y = y/(y 2 + 1) is the same as in Section IV A (with: y = A Mjv/Tjv). Even more so, numerical 


evaluations give for the N 1 -N 2 overlap parameter 5(y) the same values as in the semileptonic decays of Section IV A 
cf. Table [TV| and Figure [hf| there. 

If we assume that |-EW 2 | ~ l-^WJ (for £ = e,fi,r), i.e., Equation (92), then we have ICi ~ IC 2 = 1C, and the 
expressions for A^cp simplify significantly 


A 


(Dir.) 

7r,CP 


^ (Maj.) 
"^TT.CP 


Si „8< LNC) >lM_ =int ,|LN C ) 
(l + 6 (y) cos0( LNC )) y 

/|S eiVl | 2 sin0( LNV ) + IS^psinfl^ 0 ) 

V ISeiVrP + liW 2 


^(1 + 0(<S)) 
y 

Y-^{l + 0{8)) 
J y 


(100a) 

(100b) 


As in the case of semileptonic LNV decays of the previous Section |IV A we see that the CP asymmetry parameter 
At^cp can become appreciable and even of order one if the following two conditions are fulfilled simultaneously: (a) at 
least one of the angles 8 CC (X=LNC,LNV), defined in Equation (97), is appreciable; (b) the quantity y = AM/v/T jv is 
y ~ 1 (near degeneracy). In such cases, the estimates for the effective (true) branching ratios Br^* of Equations (63) 
and (64) would apply also to the CP-violating difference of effective branching ratios, A^ X ( lpBr^\ where (X) = 
(Dir.J^Maj.). 


V. CONCLUSIONS 

We have studied lepton number violating (LNV) semileptonic decays of charged pseudoscalar mesons, specifically 
t: ± , A , D ± , Df, B ± and , mediated by heavy neutrinos that can go on their mass shell. 

We first presented the LNV semileptonic decays of charged Kaons and of the heavier mesons D ± , Df, B ± and B±, 
in processes of the form mediated by on-shell massive neutrinos, where M is the decaying meson 

and M a correspondingly lighter meson. We estimated the branching ratios as functions of the neutrino masses and 
mixing parameters, and found the scenarios where upper limits on the mixing parameters can be obtained. We also 
studied the effect on the observability of these decays due to the long neutrino lifetime, as the secondary decay vertex 
is likely to fall outside the detector for the range of neutrino masses that are relevant to these processes. 

We then presented our corresponding study of charged pion decays, which in this case are purely leptonic since 
pions are the lightest mesons. Here we can have modes that conserve lepton number (LNC) as well as modes that 
violate lepton number (LNV), if the intermediate neutrinos are of Majorana type, while only the former modes occur 
if the intermediate neutrino is of Dirac type. However, these modes are not distinguished by the final state because 
the latter involves a standard neutrino, which is not experimentally observable. We find that it could be possible to 
discern the Majorana or Dirac nature of neutrinos if one is able to observe features in the final state distribution. 

We finally explored the possibility of observing CP violation in the lepton sector using these meson decays mediated 
by massive neutrinos on shell. The CP signal in charged meson decays is the usual asymmetry between the decays 
of opposite charge mesons. We found that leptonic CP violation may show in semileptonic LNV decays of charged 
Kaons and charged B mesons, as well as in LNC and LNV decays of charged pions, depending on the mass of the 
intermediate neutrinos. It turns out that such CP violation becomes appreciable and possibly detectable if there are 
at least two heavy neutrinos almost degenerate in mass that can go on their mass shell. The neutrino mass splitting 
that gives maximal CP asymmetries is close to the neutrino decay width. This type scenario fits well into the so 
called neutrino minimal standard model (i'MSM), which contains two almost degenerate Majorana neutrinos of mass 
near 1 GeV and another lighter neutrino of mass of order 10 1 keV, a model that can explain simultaneously neutrino 
oscillations, the dark matter and the baryon asymmetry of the Universe. 
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Appendix 


A.l. Explicit Formulas for Amplitudes of the Nj-mediated Decay AI ± —> PfufM ^ 

In this Appendix we provide, for completeness, formulas which are used in Sections [TT] and |IVAl The 
formulas were presented in Ref. .39] for the case of exchange of two different neutrinos N\ and N 2 . 
Here we present them in a slightly more general form, when the number of exchanged neutrinos is Af 
(N-\,, A(y). In Section [TI] the simpler case of Af = 1 is taken, because such a case is representative 
enough for the consideratio n o f the branching ratios there. On the other hand, the case Af = 2 (or: 
Af > 2) is taken in Section |lV| with two of the neutrinos (Ni, N 2 ) considered on-shell and almost degenerate, 
because in such a case significant CP violation effects can arise in the Majorana neutrino sector. 

The amplitude squared |T(M ± )| 2 for the decay of Figure [l] appears in the expression Equation ^ for the decay 
width T(M ± —» ifi^M : ' = ), and can be written in the form 


TV A r 

| = K 2 EE k ( j ±) kf ) *M Ni M N] 

»=i i—i 

x [P i (.D)P,(.D)*T ± (.D)T ± (.D)* + P i (C')P i (G)*r ± (C)T ± (G)* 

+ (P i (D)P J {C)*T ± (D)T±(C)* + Pi(C)Pj(D)*T±{C)T ± (D)*)] (101) 


Here 
X = 


ure 


0 


i,j = 1,...,A f are indices of contributions of the exchanges of intermediate neutrinos N i: Nj, and 
D, C denote contributions of amplitudes of the direct and crossed channels, respectively, cf. Fig- 
Further, kj^ are the heavy-light mixing factors for Nj defined in Equation (70); Pj(X ) 
(j = 1,2; A' = D,C) are the propagator functions of Nj neutrino for the D and C channel, Equation (72). K 2 
is the constant originating from the vertices and is given in Equation ([9]). These expressions appear in the normalized 
decay widths T± (XY*) in Equation |7| wh en J\f = 1, and in r±(XF*)jj in Equation (71) when J\f > 2. The quadratic 
expressions of T±(X)T±(Y)* in Equation (101) get simplified after summation over the final helicities of the leptons 
£i and £ 2 , and acquire the following form: 


T±{D)T±(D)* 

T±(C)T±(C)* 

T±(D)T±(C)* 


T±(C)T±(D)* 


8 [ M m m m'(Pi -P2)~ 2Mm(.Pi ■ PM')ip2 ■ Pm') ~ 2 Mm’(Pi • Pm)(P2 ■ Pm) 
+4(pi ■ Pm)(,P2 • Pm')(pm ■Pm ')] = T(D)T(D)* 

8[Ml I Ml I ,(p 1 -p 2 ) - 2Ml I {pi ■p M '){p2 - PM') - 2M^,(pi ■ p M ){p 2 ■ pm) 

+4(p 2 -Pm){pi ■Pm')(pm ■Pm ')] = T{C)T(C)* 

m{m'm(Pi-Pm'){P 2 ■ p M ') + Mh'(.P i • Pm){P2 -Pm) - ■ p 2 ) 

+{pm ■Pm') [~(Pi -Pm){P 2 'Pm') - (P 2 ■ Pm)(pi ■ Pm') + {Pm ■ Pm')(pi ■ P 2 )] 
T i(PM ■ Pm')c(pm,Pi,P2,Pm')} 

(T±(D)T±(Cyy = T t (D)T t (C)* = ( T T (C)T T (D)*y 


where we used the notation 


(102a) 

(102b) 


(102c) 

(102d) 


^( 91 , 92 , 93 , 94 ) = e VlTl2V3Vi (qi)rn ( 92 )^ ( 93 )^ 3 ( 94 ) 


Jr] 4 


(103) 


Here, e mmri3rii is t he to tally antisymmetric Levi-Civita tensor with the sign convention e 0123 = +1. 

The expression (102), together with the definition (71), imply for the normalized decay widths r±(Xy*)jj of 
Equation (71) various symmetry relations, namely that T±(DD*) and r±(CC'*) are self-adjoint (Af x Af) matrices, 
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and that elements of the D-C interference matrices T±(CD*) and T±(DC*) are simply related 

f {DD*)ii = (f{DD*) jt y , T(CC*)ij = (r{CC*)j 


r±(CD*)jj = (r^DC*)^' 


(104a) 

(104b) 


If the two final leptons are of the same flavor (£i = £ 2 ), one can use the property that the integration d 3 over the 
final particles is symmetric under exchange of p\ and P 2 (because = M ( 2 ), and we have the following additional 
symmetries: 


T(DD*) l3 = r (cc*)ij 

f ±(CT>*) y = T±(DC*)ij 

and the Af x Af D-C interference matrices T±{CD*) become self-adjoint, too. 
When Af = 1, as in Section [TTJ then we have in the case of i\ = l 2 

f{DD*) = f{CC*) {Af = 1 and t\ = £ 2 ) 


(105a) 

(105b) 


(106) 


A.2. Explicit Expression for the Function Q 


The expression in Equations ( |l2|) and (78) is arrived at by using in the integration over the phase space of three 
final particles [Equations ([3]) and Q], for the contribution of the N neutrino, the identity 


ds ( M(p M ) -t £i(pi)£ 2 (p 2 )M (jpm’)) = 


d 2 ( M(j>m ) -t ti{pi)Nj(p N )) dp 2 N d 2 (Nj{p N ) -t £ 2 {jp 2 )M' (p M /)) 
d 2 ( M(pm ) -t £2{p 2 )Nj{p N )) dp 2 N d 2 (Nj(pn) h(pi)M {p M ')^j 


(107a) 

(107b) 


The first identity can be used for the DD* contribution (where Pn = Pm — Pi) and the second for the CC* 
contribution (where pw = Pm — P 2 )- When one uses the identity (10) in the DD* contribution, and the analogous 
identity for the CC* contribution, the integration over dp 2 N becomes trivial, and the c^-type of integrations can be 
performed. Notice that this is equivalent to the factorization approach T{M —► IxAQBrflV —> £ 2 M ), which holds when 
N is on-shell. The obtained expression for T{DD*) is then the expression Equation ( |12| when Af = 1 [Equation (78) 
when Af > 2] with the notations (14), where the obtained function Q has the following form: 


Q{yN;ye 1 ,ye 2 ,y') = \ 1 (vn - yeJiyN - ye 2 )(i - yw - yij (1 - — + — 

2 \ VN VN 


+ [ - ye (1 + y' + 22/at - ye 1 -ye 2 )~ y\ {yw - y') + yf 2 (i - 2 in) 
+2/ti(l + vn){vn - 2 /) - yt 2 (l - vn){vn + y’)\ 


In the limit of massless charged leptons (y^ = ye 2 = 0), this reduces to 

Q{y N ; 0,0, 2 /) = \y%{l - 2/jv) f 1 - — ) 
2 V VN) 


(108) 


(109) 


A.3. Calculation of the Total Decay Width of Neutrino N 

In this Appendix, for completeness, we summarize the formulas needed for evaluation of the total decay width of a 
massive sterile neutrino N, cf. Equations {16)—(18) and Figure [2] 

In Ref. [2j (Appendix C there), the formulas for the leptonic decay and semimesonic decay widths of a sterile 
neutrino N have been obtained, for the masses Mjy < 1 GeV. For higher values of the masses Mpj, the calculation 
of the semileptonic decay widths becomes difficult because not all the resonances are known. Hence, for such masses 
the authors of Refs. [501ITT5] proposed an inclusive approach, based on duality, for the calculation of the total 
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contribution of the semileptonic decay width of TV. It consists of representing the various (pseudoscalar and vector) 
meson channels by quark-antiquark channels. This approach was applied for My > M i ss 0.958 GeV. Here we 
summarize the formulas given in Ref. 0D] for the decay width channels (cf. also: [29]). In some of these formulas, 
twice the decay width is given [2T(TV —» ...)], signalling the fact that for each possible decay of Majorana neutrino 
in charged particles, there is an equally possible decay into charge conjugate channel (something not possible if TV is 
Dirac particle). 


2 r(7v-+r/“V) 

r (tv 1//-/+) 


v ' 


\ B w\ 2 ^r^ M Nh( x e,0, x e')(l - S u >) , 

\bw\ 2 -^ m n +tu' 9 R pt) )i 2 ( 0 , , x e i) 

+ ((-?i lept) ) 2 + (^ ept) ) 2 + V(1 + 2si lept) )) h(0,x e i,x,i)] 

i 


(110a) 


(110b) 

(110c) 


In Equation (110a) factor 2 was included because for Majorana neutrino TV both decays IV —► l £ + ry and TV —> 
£ + t , ~v i > contribute (t ^ /). 

When My < M > « 0.968 GeV, the following semimesonic decays contribute, which involve presudoscalar (P) and 
vector (V) mesons: 


2 T(N ^£~P + ) 

r (tv -»■ Vi p°) 

2T(TV-> e~V + ) 
T(TV -)• veY°) 


\B m \ 2< ^M 3 N f 2 \V P \ 2 F P (xt,x P ) 

o7r 

\B m \ 2< ^M 3 N f 2 \V v \ 2 F v (x e ,x v ) 
\B tN \ 2< ^M%f 2 K 2 v {l - 4) 2 ( 1 + 24) 


(111a) 

(111b) 

(111c) 

(Hid) 


where factor 2 in the charged meson channels appears because both decays TV —)■ £ M + and TV —► £ + M contribute 
(M = P, V) if TV is Majorana. The factors V p and Vy are the CKM matrix elements involving the valence quarks 
of the mesons; and f P and fy are the corresponding decay constants, whose values are given, e.g., in Table 1 in 
Ref. [50]. The pseudoscalar mesons which contribute here are: P ± = 7 r , IV , P° = 7r°, K°, K°, 77 . The vector 
mesons which contribute are: V ± = p ± ,K* ± - V° = p°,u, K*° , K*°. If My > (=0.9578 GeV), due to duality the 

(many) semimesonic decay modes are represented by the following quark-antiquark decay modes [59i: 

2r(TV -»• i~UD) = \B eN \ 2 ^M s N \Vun\ 2 Ii(x e ,xu,x D ) (112a) 

T(N -H/ e qq) = \B m \ 2 x q , x q )+ + {9 ( r) 2 ) h(0, x q , a;,)] (112b) 


In the formulas (110)-(112) the notations xy = My /My (Y = i. vp , P. V. q) are used; and in Equation (112) we 
denoted: U = it, c; D = d, s, 6; q = u, d , c, s, b. The values of quark masses which we used in our evaluations are: 
M u = M d = 3.5 MeV; M s = 105 MeV; M c = 1.27 GeV; M b = 4.2 GeV. 

We note that in the evaluation of the total decay width Ty, the expressions (112a) and (112b) should be added 
when TV is Majorana; if TV is Dirac, the same summation should be taken, but the expressions (112a) should be 
multiplied by 1/2. The same approach is valid also in the case of summation of expressions (110) and (111). 

In Equations (110b) and (112b) there appear the following SM neutral current couplings: 


(lept) 

91 = 

+sin 2 6» w , 

Jlept) 

Sr 

= sin 2 9 W 

(113a) 

II 

cn 

1 2 . 2 

2 - 3 sm 0 W , 

to 
to ^ 

II 

2 . 2 

- - sin 9 W 

(113b) 

II 

1 1 • 2 „ 

~2 + 3 sm 

(U) 

> 9r 

= - sin 2 9 W 

(113c) 
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Further, the neutral current couplings Ky of the neutral vector mesons in Equation (11 Id) are 


Ky = — sin 2 9\y (V = p°,bj) 
o 

ny = -^ + ^sin 2 0 w (V = K*°,K*°) 

The following kinematical expressions I±, I 2 , Fp and Fy were used: 

h(x,y,z) = 12 / — (s - x 2 - y 2 )(l + z 2 - s)\ 1/2 (s,x 2 ,y 2 )\ 1/2 (l,s, z 2 ) 

J(x+y ) 2 s 

h{x,y,z) = 24 yz [ -(l + x 2 -s)\ 1/2 (s,y 2 ,z 2 )\ 1/2 (l,s,x 2 ) 

J(y+z ) 2 s 

Fp(x, y) = A 1 / 2 (l, x 2 , y 2 ) [(1 + x 2 )(\ + x 2 — y 2 ) — 4x 2 ] 

Fy(x, y) = A 1 / 2 (l, x 2 , y 2 ) [(1 — x 2 ) 2 + (1 + x 2 )y 2 — 2y 4 ] 


(114a) 
(114b) 

(115a) 

(115b) 

(115c) 
(115d) 

with A function defined in Equation ( |46a[ ) . These formulas allow us to obtain the total decay width T(N —>• all) as a 
function of Mpf. Using these formulas, we evaluated the coefficients Afm, appearing in Equation ( p~8|) at the mixing 
terms \Bin\ 2 , and presented them in Figure [ 2 ] as a function of Mjv for the cases of Majorana and Dirac neutrino N. 
One may notice a small kink in the curves of Figure [ 2 ] at Mjv = M^ (=0.9578 GeV). This kink appears because 
at Mpj > M 1 the use of duality is made (the replacement of the semileptonic decay channel contributions by the 
quark-antiquark channel contributions). We can see that the duality works quite well at Mjy > , with the possible 

exception for the case t = r because r lepton has a large mass. 


A. 4 . Explicit Amplitudes for Nj-Mediated Decays 




In this Appendix we summarize, for completeness, formulas needed in Sections |III| and |IV B[ These formulas were 
derived and presented in Ref. [381 , for the case of exchange of two different neutrinos and N 2 ■ Here we summarize 
them in a slightly more general form, for the case of A f different neutrinos Nj (j = 1,..., A f). In Section 


III 


the 


simpler case N = 1 is taken, as it is sufficiently representative for the branching ratios considered there. In Section IV 


the case JV = 2 (or: J\f > 2) is considered, with two (on-shell) neutrinos N\ and N 2 almost degenerate, as in such a 
case significant CP violation effects can occur in the neutrino sector. 

The squared amplitude |7^^| 2 for the Nj -mediated leptonic decays of neutrinos, appearing, for example, in Equa¬ 
tion ( [36| ) (where X = LNC, LNV), is a combination of contributions from the two channels D (direct) and C (crossed) 
(c/. Figures [8] and [9]), and, in general, of the contributions of Af neutrinos Nj 


N a r 


rt,?i 2 = *2££*2 ) '*£ ) 


*=1 0= 1 


P^ X) (D)Pj X ^ (D)*T^ x 2 (DD*) + Pl x) (C)Pj x) (C)*t( x 2(cc*) 

^ X) {D)pj x) (C)*T^ X 2(DC*) + Pj; X \c)P^ X \D)*Ti x 2{CD*) 


(116) 


The constant K 2 is given in Equation (42), and the mixing factors k^2 are 


, (LNV) _ r2 
K j,+ ~ & 


eNi 


7 (LNC) _ R R * 
k j,+ -FeNjBfjNj 




k'Cf = 


(117) 


In the case of J\f = 1, these coefficients are in Equation (39). P> (Y) are the A^-propagator functions [when 
J\f = 1: Wpropagator functions of Equation ©1 of the direct and crossed channels (Y = D,C). 

Explicit expressions for the direct ( DD *), crossed (CC*) and direct-crossed interference (DC* and CD*) terms 


[t22(F>D*), t22_(CC*), t2%(DC*), tI X 2(CD *)\, appearing in Equations (116), get simplified when summed over 


the helicities of all the final leptons. In the case of the X = LNV processes (c/. 


’igure [9]) they acquire the following 
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form: 


r£f v) (£>£>*) = 256(p2 • p„) [-M 2 (p x • p„) + 2(p! • p ff )(p„ • p w )] = 4 LNV )(£D*) , 

T^ V) (CC*) = 256(p! • p„) [-Af 2 (p 2 • p„) + 2(p 2 • p„){ Vfl ■ p n )\ = ri hNV \CC*) , 

Ti™ V) (DC*) = 12s|(pi • p„) [AT 2 (p 2 ■ p„) - 2(p 2 • p ff )(p M • p*)] 

+(P2 • Pt/) [M 2 (pi • p M ) - 2(pi • p*.)(p M • p*-)] 

-(Pi • P 2 ) [Af 2 (p„ • p M ) - 2(p„ • p 7r )(p /i • pj] | 

T*j - (pi * P7r)^(P2 1 Pisi Pfi 5 P-7T) + (P2 ' PttMPI, P^, P M , Ptt) 

-(Pi/ •P7r)e(Pl,P2,P jU ,Pir) - (Pm •P7r)e(Pl)P2,Pi/,Pir) 
r (LN v) (c ^) = (tCLNV)^)^ =T (LNV) ( ^ } 


where we used the notation Equation (|103j) for e. 

In the case of X = LNC processes (c/. Figure [8]), the expressions are 


T^ C \DD*) = 256(p /i • p v ) 


(pi • P 2 ) (M* - Af 2 AL 2 - 4M 2 (pr • p ff ) + 4(pi • p w ) 2 


+2M e 2 (p 2 • Ptt)(M 2 - pi • p w ) 


= T^ LNC) (£)£>*) 


T^ c )( C 'C*) = 256(p At • p„) 


(Pi • P 2 ) (Af* - Af 2 M 2 - 4A/ 2 (p 2 • p ff ) + 4(p 2 • p w ) 2 


+2AF e (pi • p 7r )(M 7r - p 2 • p,r) 


= r^ LNC )((7C*) 


t| L ± C) PC*) = 256 (p M -p„) 


(pi • p 2 )(Af 2 - 2pi • p w )(Af 2 - 2p 2 • Pvr) 


(118a) 

(118b) 


(118c) 

(118d) 


(119a) 


(119b) 


(119c) 


+M 2 (—2(pi • p T ) 2 - 2(p 2 • p w ) 2 + Af 2 (pi + p 2 ) • p w + AP 2 M e 2 ) 

= T^ lng) (Z)C*) , 

T( lnc )(CD*) = ( T^ nc \DC*)Y (119d) 

These expressions appear in the definition of the normalized (i.e., without explicit mixing dependence) decay width 
matrices f { *l(YZ*) i:j (X = LNV, LNC; Y, Z = D, C; i,j = 1,..., TV) 


-(X) 

■ 7 r,± 


(YZ*)ij = K 2 


1 1 


2! 2A4 (2tt) 8 


d 4 P> ’{Y)P^ K) (Z)* T±'(YZ*) 


( 120 ) 


where -Pj X ' 1 are the propagator functions of neutrino Nj (with mass ). cf. the definitions (|4l|) written when TV = 1. 

When TV = 1, the definition (120) reduces to the definition (|40| in Section [ill A 


When we use the symmetry of the d 4 integration under the exchange p 4 O p 2 (we note: 
Mi = M 2 = M e in our considered case), this leads to the following identities: 


rW(£>D*)y = Y^XCC*)^ , r ( ^\DD*) ji = (Ti x \DD*) ij y (121a) 

rg(C£>*)« = f i X l(DC*) l0 = (f^ (<?£*),,)* ( 121 b) 

In the case of TV = 1 this reduces simply to 

f W(DD*) = f l x) (CC*) (122) 
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-wv\ ~/va 

A.5. Explicit Expression for Tn ; and dYi /dE M 


for 7r =t: —> e ± e ± pf f u with On-Shell N 


Equation (451 refers to the expression obtained by performing the integration (36) over the phase space of the four 
fina l particles [cf. Equation (37)], of the integrand written explicitly in Appendix 4. In the integration, the on-shellness 
(43) is assumed, which makes the integration over p 2 N trivial. At the final stage of integration, the differential decay 
width dT ( x )/dEp over the muon energy E^, in the rest frame of the N neutrino, is performed. The expressions for 
df^/dE^ were written in Refs. E31SH1, and we write them down here for completeness. In the case of X = LNV it 


is 


C?r(LNV) e ± e ±fjT v ) 


dE,. 


Ki 


2(27t) 4 Y n M% 


A Ml) x [MlMl - M% + M e 2 (M 2 + 2 M 2 N - M 2 )] 

(M 2 N + Ml-Mp \ 


xE ~ 2MnE » + M l ~ M « ) 2 ( 

" (M 2 -2M N E ll +MZ) ^ 


< E^ < 


2 M„ 




(123) 


The integration of this expression over E M can be performed explicitly (in Ref. m it was performed only in the 
limit M e = 0). The result is Equation (45) with notations (46), where the function E(xp, x e ) was obtained in Ref. [38j . 
We write it down here again, for completeness. 

Eix^Xe) = | A 1,/2 (l, X/j,, x e ) [(1 + - 8Xfj, + a; 2 ) - x e (7 - 12x M + 7x 2 ) 

-7x 2 e (l + x^) + x 3 e ] — 24(1 — a: 2 )a; 2 In2 


+12 


- x 2 (l - x 2 e )\nx^ + (2x 2 - x 2 (l + x 2 ))ln(l + 0 ,^ + A 1/2 (l, x^, x e ) - x e ) 


„ 2n „ 2 m„ f (l ~ x^) 2 + {l - x fl )\ 1/2 (l,x fl ,x e ) - x e (l + x^) 


+X e (l — x u) l n 


(124) 


We can obtain the LNV canonical differential decay width, according to Equation (54), from the normalized 


differential decay width Equation (123). For this, it turns out to be convenient to use the following identity: 

Ki 


2K 


= 96 


1 


T(tt± all) (2 7 r) 4 r JV M3 - M 2 ) 2 (1 + Sg v ) 


(125) 


which is obtained by using Equations (42), (16) and (H7k and (50). This then gives us 


dBr 


(LNV) 


dE a 


= 48 


= 2 


k. dr( LNV ^(7r d 


o± p± ,,=F 




Yr. 


dE a 


-M 2 ) 2 ( 1 + Sg*) M ^ X ~ M n + + * M n ^e)} 


-{M 2 n - 2M N E fl + Mi - Ml f 


* E »\' E i - M r N {M r^ NE ;: M ^ > ^ 


(126) 


where 


(^) 


max — 


(M 2 N + Ml~Mj) 

2Mjv 


(127) 
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The LNC canonical differential decay width turns out to be 

K df^ G \^ ^ e ± e ± pAv) 


d^ C) 


dE„ 


= 2 - 


r„ 


dE u 


-A 1 / 2 (M 2 , Ml, Ml) - 


+ W -[Ml 

x js y/(E*-M*)M N [(2 Eu. - M n )M n - M 2 + M 2 ] 2 
x [M 2 M 2 - M% + Ml (Ml + 2M 2 N ) - Ml] 


■M N (-2E^ + M N )Y 


%ElM 2 N - 2M 2 Mjv(M 2 + + 2M 2 ) - 2£ 2 M W (5(M 2 + M 2 ,) + M, 


+S M (3M^ + 10M 2 M 2 + 3M* + 3M 2 (M 2 + M 2 ,)) 


(M m < < (Efi ) m ax) 


(128) 


It turns out that, upon integration of this expression over E M , we obtain the same result as in the X = LNV case, i.e., 
Equations (53c) with (46) and ( |124 ), or equivalently, Equations (45) with (46) and (124). We must add that we found 
a typographical error in Equation (A.16) of Ref. i[37], where E 2 must be replaced by 2 E\, and in Equations (B.lc) 
and (B3) of Ref. [35|, where (E ^) max should read (Mj^ + Ml — Ml)/(2M N ). 

In the limit M e = 0 (which is a good approximation), the canonical differential decay widths (126) and (128) get 
simplified 


d&i LNV) 


dEp 

dB?i LNC) 


dE u 


M e —0 


M e = 0 


- 2M » E " + "A 


48(M 2 - M%) 2 /^2~M2 

M%M*(M*-M*)*V " " 


-En(Mfj + Ml) - -Mtv(2E 2 + M 2 ) 


(129a) 

(129b) 


The full (integrated) canonical branching ratio in the M e = 0 limit is obtained by taking the x e = 0 limit of 
Equation ([53c| 


Br„ 


1 


lM e =0 


M 2 n (Ml - M*f ( Ml 

2 Ml (Ml-MlY J \Ml 


(130) 


where the function / is written in Equation (481. 

A.6. Delta Function Approximation for the Imaginary Part of the Propagator Product 


In this Appendix we investigate the expression for the imaginary part of the propagator product, lm(Pi(D)P 2 (D )*), 
Equation (79a I of Section IV A For convenience we introduce in this Appendix the following simplified notations x, 
M 2 , A and £: 


x = p N , 


M 2 = Ml x 


A = A Ml = Ml 2 — Ml x 


■ Ni 


= £Tjv, Fn 2 = (2 - O^N 


(131a) 

(131b) 

(131c) 


We note that A > 0 by convention; and 0 < £ < 2. Further, T,y, + Pn 2 = 2T n, in accordance with the definition 
of Tat Equation (82). Since we always have Tjv. -C Mw (the neutrinos Nj are sterile), the relation (80) holds, i.e., 


r N m n 


(X-MIV+TIM* 


= 7 T§(x-Ml) 


(132) 


We can write the right-hand side of Equation (79a) for lm(Pi(D)P 2 (D)*) as 


Im(P 1 (D)P 2 (D)*) = n 1 +n2 


(133) 
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where 1Z\ and 77-2 can be written, in our notation, as 

(x - M 1 2 3 4 5 6 7 8 9 10 11 ){2 - O^nVM 2 + A 


n, = 


[{x - M 2 ) 2 + £ 2 r 2 n M 2 } [{x - M 2 - A) 2 + (2 - 0 2 r^r(M 2 + A)] 


= 771 x —<J(a; ~ M 2 ~ A) 


n 2 = - 


£T n M(x - M 2 - A) 


[(* - M 2 ) 2 + £ 2 r^M 2 ] [(a: - M 2 - A) 2 + (2 - ^) 2 r 2 r (M 2 + A)] 


= V 2 x ~^6{x - M 2 ) 


(134a) 

(134b) 

(134c) 

(134d) 


In Equations (134b) and (134d), the identity (132) was used, and we introduced two (dimensionless) parameters r/j 
(j = 1,2). We want to obtain these two parameters r)j. They can be obtained by integrating analytically the explicit 
expressions (134a) and (134c) for lZj(x) over x. For example, integration of IZi(x) gives 


p+oo 


dx 


(x - m 2 )(2 - or n VM 2 + A 


7T A 


[(a; - M 2 ) 2 + Or 2 N M 2 ] [(x - M 2 - A) 2 + (2 - £) 2 T 2 N {M 2 + A)] (A 2 + 4T 2 N M 2 ) 


where 


Mi = -M 2 


= M 


(2 - ^(2 - 0) + £(2 - OVl + A/M 2 ] + -(2 - C) 2 A 
A / A 2 \ 

! +(1-^/2)-^2 +0 ( 


M 4 ) 


(135) 

(136a) 

(136b) 


Therefore, in the c ase of near degeneracy (A <C M 2 ) we have M 2 — M 2 . If we now use in the integration o ver dx 
the expression (134b) instead, take into account M 2 = M 2 in the case of near degeneracy, and compare with (135), 
we obtain the following expression for the parameter rji by comparison with (135): 

A 


1 

?/l A ~ (A 2 + 4T 2 n M 2 ) 


y 


= 2 

y 2 


1 


(A < M 2 ) 

A 


y = 


(2 MT n ) 


, A < M‘ 


(137a) 

(137b) 


where in Equation (137b) we use the usual notation in this paper y = (Mn 2 — MjVjJ/Tjv = A/(2 MT^). Here we note 
that A = (M^ 2 — MfjO = (M N2 — M Nl )2M Nl in the case of near degeneracy A M 2 = 

Doing the same procedure with the quantity 1Z 2 , we obtain for rj 2 the very same result as for 771 


y 2 + 1 


(A < M 2 ) 


(138) 
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